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ABSTRACT 


A  model  is  developed  for  relating  continuous  review 
inventory  policies  for  repair  parts  to  system  availa¬ 
bility.  The  system  consists  of  S  identical  unit  systems, 
each  of  which  is  a  series  of  k-out-of-n  structures.  Unit 
system  states  are  zero  or  one.  An  optimal  cannibalization 
policy  is  assumed.  Under  this  assumption  the  number  of  unit 
systems  up  is  always  the  maximum  possible  for  any  given 
vector  of  backorders  for  the  N  part  types  in  the  system. 
The  distribution  of  backorders-  under  a  (Q,r)  policy  with 
Poisson  demands  for  each. -part  type  is  used  to  derive 
expressions'  for  system  availability  as  functions  of  the 
Q,  r  vectors .  For  simplicity  it  is  assumed  that  order 
quaj.ti-ties  *re  set  by  an  ope'rating  level  in  terms  of  days 
of  supply/.  A  numerical  technique  is  presented  for  finding 
the  vector  of  reorder  points  (safety  levels)  which  mini¬ 
mizes  the  expected  cost  of  on  hand  inventory  subject  to  one 
of  two  alternative  availability  constraints. 
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I .  INTRODUCTION 


A.  BACKGROUND 

In  the  management  of  logistic  support  for  weapons 
systems  a  significant  problem  in  resource  allocation  is 
the  determination  of  optimal  inventory  policies  for  repair 
parts.  Inadequate  stocks  of  repair  parts  result  in  low 
system  availability.  Stocks  sufficient  to  insure  with 
high  probability  that  any  part  will  be  immediately  avail¬ 
able  when  needed  my  tie-up  more  resources  than  are 
justified. 

An  indication  of  the  concern  of  the  logistics  manager 
with  this  problem  is  the  following  statement  from  the 
revised  General  Objective  No.  1  for  the  Navy  Supply  Systems 
Command:  "For  technical  material,  optimum  support  is  that 

which  minimizes  downtime  of  weapons  systems  due  to  lack  of 
repair  parts  and  components."  [Ref.  1]  In  order  to  accom¬ 
plish  this  objective,  techniques  are  needed  for  relating 
repair  parts  inventory  policies  to  weapons  systems 
availability . 

The  mathematical  models  designed  to  deal  with  the  impact 
of  repair  parts  shortages  have  usually  been  formulated  as 
a  minimization  of  total  cost,  where  the  total  cost  includes 
ordering,  holding,  and  backorder  costs.  In  theory  the 
backorder  costs  are  a  measure  of  the  impact  of  shortages  on 
system  availability.  A  difficulty  in  practice  is  that 
backorder  costs  are  hard  to  estimate.  One  technique  for 
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finding  a  value  of  the  shortage  cost  is  to  impute  it  based 
on  a  required  maximum  probability  pout  of  being  out  of 
stock.-  See  Deemer  and  Hoekstra  [Ref.  2,  p.  5]  for  a  dis¬ 
cussion  of  this  technique.  The  required  value  of  PQut  is 
determined  by  management  judgement.  This  judgement  pre¬ 
sumably  includes  some  intuitive  consideration  of  the  effect 
of  repair  parts  shortages  on  system  availability. 

The  total  cost  minimization  formulation  does  not,  how¬ 
ever,  permit  any  direct  correlation  of  inventory  policies 
and  system  availability.  Inventory  policies  computed  using 
that  formulation  are  determined  for  each  part  independent 
of  all  other  parts.  System  availability  is,  however,  a 
function  of  (Q,r)  =  (Qx »rx) , . . . , (QN ,rN) .  Therefore,  if 
we  wish  to  correlate  inventory  policies  to  system  avail¬ 
ability,  we  need  a  multi-item  model  rather  than  one  which 
deals  with  each  item  one  at  a  time. 

There  do  exist  a- number  of  techniques  for  explicitly 
correlating  system-  availability  and  single  period  inventory 
policies.  One  of  these  is  the  optimal  redundancy  approach 
discussed  in  Chapter  6  of  Barlow  and  Prochan  [Ref.  3] .  The 
problem  discussed  there  is  one  where  at  the  beginning  of 
the  period  there  is  a  quantity  n..  of  each  part  type  j 
on  hand,  and  resupply  is  not  available  until  the  end  of 
the  period.  The  object  is  to  determine  the  value  of  each 
nj  so  that  system  availability  is  maximized  subject  to 
some  cost  constraint  or,  alternatively,  to  minimize  the 
cost  of  achieving  a  required  availability.  The  cost  of  the 
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parts  inventory  may  be  interpreted  as  dollars,  weight, 
volume  or  some  other  measure  of  the  amount  of  a  ’-'resource" 
which  is  used  up  by  the  inventory. 

B.  GENERAL  NATURE  OF  THE  MODEL 

The  approach  taken  here  is  to  apply  the  redundancy 
optimization  ideas  to  a  continuous  review  inventory  sit¬ 
uation.  In  particular  the  model  addresses  the  problem 
of  optimal  parts  inventory  policies  in  support  of  a  weapons 
system  at  the  direct  support  echelon.  The  objective 
function  to  be  minimized  is  the  expected  cost  of  on  hand 
inventory.  The  constraint  is  a  system  availability  re¬ 
quirement,  which  may  be  in  one  of  two  alternative  forms. 

The  first  is  that  the  expected  number  of  unit  systems  up 
must  be  greater  than  or  equal  to  a  required  fraction  of 
the  total  number  of  unit  systems.  The  second  form  of 
the  constraint  is  that  the  probability  that  at  least  k 
unit  systems  are  up  must  be  equal  to  or  greater  than  a 
required  assurance  'level. 

The  model  is  an  idealization  of  the  repair  parts 
supply  support  furnished  by  the  direct  support  maintenance 
unit  for  the  population  of  a  particular  major  item  in  an 
operational  Army  unit.  The  system  consists  of  S  identical 
unit  systems,  which  might  be  tanks,  aircraft,  howitzers, 
or  some  other  set  of  equipments  of  the  same  make  and 


model. 


C.  .ASSUMPTIONS  OP  THE  MODEL  AND  THEIR  IMPLICATIONS 

Before  turning  to  the  formal  description  of  the  model 
we  shall  discuss  some  of  the  important  assumptions  and 
their  implications  here. 

1.  The  Inventory  System  is  Single  Echelon 

To  simplify  the  problem  we  assume  that  the  supported 
military  unit  does  not  carry  a  stock  of  repair  parts,  but 
that  parts  are  immediately  available  from  the  direct  support 
maintenance  unit  if  there  is  stock  on  hand  at  the  time  a 
demand  occurs.  The  inventory  system  thus  is  single  echelon. 

2 .  Procurement  Lead  Time  is  Constant 

The  procurement  lead  time  (order  and  ship  time)  is 
assumed  to  be  constant  for  each  part  type where  part  type 
indicates  a, ^particular  Federal  Stock  Number.  Lead  times  may, 
however,  be  different  for  different  part  types.  We  further 
assume  that  the  source  of  supply  is  never  out  of  stock. 

Thus,  the  model  ignores  variability  of  lead  times  and  the 
possibility  of  reducing  lead  times  by  using  a  higher  requi¬ 
sition  priority. 

3 .  Unit  System  State  Values  are  Zero  or  One 

The  supported  system  consists  of  S  identical  unit 
systems.  Each  unit  system  is  assumed  to  exist  in  one  of 
exactly  two  states  depending  upon  the  states  of  the  in¬ 
stalled  parts:  either  it  is  up,  i .e. ,  it  is  fully  capable 
of  operating  satisfactorily;  or  it  is  wn-,  i.e.,  totally 
ineffective.  This  assumption  means  tha**  'the  possibility 
that  the  unit  system  may  be  partially  effective  is  not 


14 


considered.  Actual  equipment  does  exist  _n  partially 
effective  states.  For  example,  if  the  50-caliber  machine 
gun  on  a  tank  is  not  working,  the  capability  of  the  tank 
is  reduced  in  some  of  its  roles,  but  it  is  obviously  not 
totally  ineffective  if  all  other  components  of  the  tank 
are  working.  The  zero-one  assumption  is  made  for  mathemati¬ 
cal  simplicity.  The  two  state  categorization  is  not  too 
different  from  the  equipment  serviceability  code  (ESC) 
ratings  currently  us,ed  within  the  Army  for  determining  the 
materiel  readiness  of  Army  units.  Possible  ESC  ratings  are 
green,  meaning  fully  operational  and  capable  of  operating 
in  combat  for  60  days;  amber,  meaning  fully  operational  and 
capable  of  operating  in  combat  for  30  days;  and  red,  mean¬ 
ing  not  fully  operational.  We  have  simply  dropped  the  amber 

* 

category . 

4 .  Component  Part  State  Values  are  Zero  or  One 
Similar  to  the  zero-one  assumption  for  each  unit 

system  two  possible  states  are  assumed  for  each  installed 
part:  either  it  is  working  or  it  has  failed.  Thus  the 
part  is  100  per  cent  effective  or  totally  ineffective. 

5 .  The  Unit  System  is  a  Series  of  k-out-of-n  Structures 
As  mentioned  in  Section  3  above,  the  state  of  a 

unit  system  is  determined  by  the  states  of  its  component 
parts.  If  there  are  a.,  parts  of.  type  j  initially 
installed  on  each  unit  system,  it  is  assumed  that  there  is 
a  number  bj  equal  to  or  less  than  a^  such  that  at  least 
bj  parts  of  type  j  must  be  working  if  the  unit  system  is 
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to  be  up.  If  the  number  of  working  type  j  parts  is  less 
than  bj,  the  unit  system  will  be  down.  Further  we  assume 
that,  if  the  number  of  working  type  j  parts  is  equal  to 
or  greater  than  bj  for  all  j ,  then  the  unit  system  will 
be  up.  These  assumptions  may  be  summarized  by  saying  that 
the  unit  system  is  a  series  of  k-out-of-n  structures,  as 
defined  by  Birnbaum,  Esary  and  Saunders  [Ref.  4,  p.  58]. 

Our  b j  corresponds  to  these  authors '  k  and  our  aj  to 
their  n. 

6 .  The  Number  of  Unit  Systems  Up  is  the  Minimum  of  the 

Number  Up  with  Respect  to  Each  Part  Type 

If  2  is  the  number  of  unit  systems  up  at  some 
arbitrary  time  and  Zj  is  the  number  which  would  be  up  if 
all  non-working  parts  other  than-  type  j  were  replaced 
by  working  parts,  then  Z  is  called  the  state  of  the 
system,  and  Zj  is  called  the  state  of  the  system  with 
respect  to  part  type  j .  We  assume  that  at  any  time  t  the 
state  of  the  system  Z  equals  the  minimum  over  j  of  the 


7 .  An  Optimal  Cannibalization  Policy  is  Observed 

Perhaps  the  most  significant  assumption  of  the  model 
c.yeloped  here  is  that  a  policy  of  optimal  cannibalization 
is  observed.  By  optimal  we  mean  that  given  any  vector 
W  =  (W  , ...,WM),  where  W.  denotes  the  number  of  working 

.IN  3 

parts  of  type  j ,  after  the  cannibalization  operation  the 
maximum  possible  number  of  unit  systems  will  be  up.  The 
mathematical  structure  of  the  supported  system  is  an 
adaptation  of  the  -structure  of  the  cannibalization  model  of 
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Hirsch,  Meisner  and  Boll  [Ref.  5].  The  reason  for  using 
this  structure  is  that,  given  any  vector  Y  =  (Y^ , . . . ,  Y^) , 
where  Yj  denotes  the  number  of  backorders  of  part  type 
j,  there  is  a  unique  number  systems  up.  This  would  not  be 
the  case  if  cannibalization  were  not  allowed.  Consider, 
for  example,  a  system  in  which  each  unit  system  has  two 
parts  of  type  j  installed,  and  both  of  these  must  be 
working  if  the  unit  system  is  to  be  up.  Suppose  there  are 
two  backorders  for  this  part  type.  This  means  that  there 
are  two  non-working  parts  of  type  j  in  the  system.  If 
one  of  these  were  on  one  unit  system  and  the  other  were 
on  a  different  unit  system,  then  the  number  of  unit  systems 
down  for  this  part  type  would  be  two.  If,  however,  both 
were  on  one  unit  system,  then  only  one  unit  system  would 
be  down  for  this  part  type.  Under  a  policy  of  optimal 
cannibalization  there  would  always  be  only  one  unit  system 
down  for  this  backorder  situation.  Since  the  probabilities 
for  the  number  of  systems  up  used  in  this  model  are  com¬ 
puted  based  on  the  probability  distributions  of  the  Y  ^ , 
the  cannibalization  assumption  is  an  essential  feature  of 
the  model.  The  degree  to  which  cannibalization  is  actually 
used  to  increase  system  availability  in  practice  depends 
upon  command  policy  and  the  practicality  of  taking  parts 
from  one  unit  system  to  make  another  operational.  If  the 
supported  unit  is  spread-out  geographically,  it  may  not  be 
feasible  to  cannibalize.  Also  the  maintenance  effort  in 
removing  and  replacing  parts  under  the  cannibalization 
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operation  may  overload  the  capability  of  the  maintenance 
organization.  Nevertheless,  some  cannibalization  may  be 
the  best  strategy  for  optimizing  system  availability  when 
there  are  constraints  on  the  quantities  of  repair  parts 
the  .maintenance  unit  can  carry. 

8.  The  Demand  Distribution  for  Each  Part  Type  is  Poisson 
The  demand  distribution  for  each  part  type  is  assumed 
to  be  Poisson  with  mean  equal  to  the  product  of  the  failure 
rate  for  the  part  type  times  the  number  of  parts  of  that 
type  in  the  system.  The  failure  rate  'of  each  part  type  is 
assumed  to  be  the  same  in  all  of  its  applications.  The 
Poisson  assumption  implies  that  the  mean  time  between 
failures  of  a  given  part  type  is  exponentially  distributed. 

It.  is  true  that  the  exponential  distribution  may  be  a  poor 
fit  for  the  mean  time  between  failures  where  an  individual 
part  is  installed  and  replaced  with  a  new  part  of  the  same 
type  immediately  upon  failure,  particularly  when  the  part 
is  subject  to  wearout,  Cox  [Ref.  6,  p.  77]  indicates  how¬ 
ever,  that  the  pooled  output  of  a  number  of  renewal  processes 
tends  to.  have  the  properties  of  a  Poisson  process  as  the 
number  of  renewal  processes  being  pooled  gets  large.  Con¬ 
sequently,  the  Poisson  assumption  may  be  a  good  approximation 
when  the  total  number  of  parts  of  a  given  type  in  the 
system  is  large.  Another  implication  of  the  assumption  is 
that  the  distribution  of  demands  is  no,t  affected  by  the 
number  of  unit  systems  down.  This  in  turn  implies  that  an, 
inoperable  unit <system  continues  to  generate  part  failures 
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even  after  it  has  gone  down.  Errors  introduced  by  this 
assumption  will  not  be  too  great  if  the  availability  rate 
is  reasonably  high  or  if  failure  rates  are  estimated  on 
demand  histories  over  a  period  when  the  availability  rate 
was  about  the  same  as  the  .required  availability  rate  for 
the  period  for  which  inventory  policies  are  being  computed. 

9 .  Demand  Rates  are  Linear  with  Activity  Levels 

The  parameter  in  the  demand  distribution  is  assumed 
to  be  a  time  rate.  As  discussed  by  Soland  ,[Ref.  7,  p.  45], 
the  natural  parameter  of  the  demand  process  may  be  miles 
driven,  rounds  fired,  or  some  other  measure  of  the  usage 
or  activity  level  of  the  system.  The  device  used  in  the 
model  to  account  for  this  fact  is  an  "activity  level" 
multiplier  La  for  the  pooled  demand  rates  for  each  part 
type.  Suppose,  for  example  that  rounds  fired  were  the  nat¬ 
ural  parameter  for  the  demand  process.  Further  suppose 
that  a  certain  number  of  rounds  were  fired  on  the  average 
by  each  unit  system  per  month  during  the  period  when 
demand  data  were  accumulated  for  estimating  failure  rates. 
Now  suppose  that  during  the  period  for  which  inventory 
policies  are  being  computed  that  the  programmed  number  of 
rounds  per  unit  system  per  month  is  doubled.  Then  the  L_ 

Q 

would  be  two,  and  the  failure  rate  estimates  for  each  part 
type  would  be  doubled.  Use  of  this  device  assumes  that 
failure  rates  are  linear  with  the  level  of  usage.  This 
assumption  is  not  in  general  true,  but  any  more  realistic 
means  of  handling  this  problem  would  complicate  the  model. 
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For  part  types  net  subject  to  aging,  this  assumption  is 
probably  fairly  good. 

10 .  Parts  are  Replaced  Immediately  upon  Failure 

This  model  assumes  that  each  part  is  replaced 
immediately  upon  failure  with  zero  replacement  time  if 
a  spare  part  is  on  hand  at  the  time  of  failure.  Also,  if 
backorders  exist  for  a  particular  part  type,  the  failed 
part  will  be  installed  immediately  upon  receipt  by  the 
maintenance  unit.  Of  course,  these  assumptions  are  not 
true  in  the  "real  world."  To  account  for  the  fact  that 
repair  times  are  finite  the  following  procedure  could  be 
used.  First  estimate  the  availability  which  we  would 
expect  to  achieve  if  there  never  ware  any  shortage  of  repair 
parts.  Call  this  estimate  the  availability  with  respect 
to  maintenance.  Call  the  estimated  availability  predicted 
by  this  model  the  availability  with  respect  to  supply 
support.  Then  by  taking  the  product  of  these  two  estimates 
we  should  have  a  fair  estimate  of  the  overall  availability 
to  be  expected  given  the  inventory  policies  and  the  repair 
capabilities  of  the  maintenance  unit.  In  any  case,  the 
reader  should  keep  in  mind  that  availability  values  computed 
with  this  model  need  to  be  corrected  for  down  time  due  to 
non- zero  repair  times. 
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II.  MATHEMATICAL  MODEL 


A.  FORMAL  DESCRIPTION  OF  THE  SYSTEM  STRUCTURE 

The  system  consists  of  S  identical  unit  systems.  On 
each  unit  system  k  there  are  n  loci  where  parts  are 
installed.  Let  the  denote  the  ith  locus  on  the  kth 

unit  system  and  let  the  set 


A* =  u*iCi 


denote  the  set  of  loci  on  the  kth  unit  system.  Let  the  set 


A  °  '-J  Ak 

k=l  K 


denote  the  set  of  all  loci  on  the  S  unit  systems.  Let 
X^  be  a  random  variable  indicating  the  state  of  the  ith 
locus  on  the  kth  unit  system,  such  that 


X,  .  =  1  if  locus  Afci  contains  a  working  part  of 
the  type  required  in  that  locus, 

=  0  if  locus  Afc-[  fails  to  contain  a  working 
part  of  the  required  type. 


How  Xk  =  (Xkl . Xkn  )  is  a  vector  of  zeros  and  ones  which 

describes  the  state,  of  the  kth  unit  system.  The  possible 
values  of  X^  correspond  to  the  vertices  of  the  unit  cube 
in  Euclidan  n-space.  Let  denote  the  set  of  all  pos¬ 

sible  values  of  X^.  The  size  of  ^  is  2n.  Let  the 
random  vector 


X  =  (x1,,. . ,xs) 
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be  an  nS-component  vector  indicating  the  state  of  the 
entire  system.  Let  %  denote  the  set  of  all  possible 

ri  C 

values  of  X.  The  size  of  ^  is  2n  .  Let  the  set 

n 

$  =  {FSN  . } 

J  j=l 


represent  the  set  of  different  part  types  installed  in  the 
n  loci  of  each  unit  system.  -Further  let  G(FSNj)  be  the- 
subset  of  which  has  as  elements  all  those  loci  require 

ing  part  type  FSN..  It  is  assumed  that  if  locus  A,  .  eG(FSN .), 

1  Kl.  J 

then  X,.  =1  if  and  only  if  a  working  part  of  type  FSN. 

Kl  3 

is  installed  in  that  locus.  This  assumption  implies  that 
different  FSN's  are  not  substitutes  for  one  another.  Let 
the  size  of  G  (FSN ^ )  =  a  ^  ,  i.e.,  the  number  of  applications 

of  FSN j  on  each  unit  system  is  a^.  Assume  that  for  each 
j  there  exists  a  number  b^  such  that 


1.  1  <  bj  <  aj; 

2.  If  the  number  of  working  parts  of  type  FSN ^  in¬ 
stalled  on  unit  system  k  is  equal  to  or  greater 
than  bj,  the  unit  system  will  not  fail  due  to 
FSN  j ; 

3.  If  the  number  of  working  parts  of  type  FSN^  in¬ 
stalled  on  unit  system  k  is  less,  than  b^,  the 
unit  system  will  fail  due  to  FSN ^ ;  and 

4.  b.  is  independent  of  b.  for  all  i  ^  j. 

D  1 

Let  the  state  of  the  unit  system  with  respect  to  FSN-j , 


denoted 


be  defined  as  the  state  in  which  all  loci 
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not  in  the  set  G{FSNj)  have  working  parts  installed. 

In  other  words  U,  .  describes  the  state  of  the  unit 
system  if  FSN  j  were  only  part  type  subject  to  failure. 
If  is  a  random  variable  indicating  the  number  of 

working  parts  of  type  FSNj  installed  on  unit  system  k, 
then 


Ukj  *  1  if  Wkj  ibj' 
"  0  otherwise. 


Now  let  the  state  of  the  system  with  respect  to  FSN.. ,  • 

S  3  s 

denoted  Z . ,  be  defined  as  £  U,  . .  Z  .  indicates  the  num-» 

1  k=i  K3  3 

ber  of  unit  systems  which  would  be  operational  for  a  given 
set  of  values  of  W^  j  ,  k  =  1,...,S  if  the  only  fallable 
part  type  were  FSN  ^ . 

For  a  given  value  of  the  system  state  vector  X  it 
may  be  possible  to  increase  the  value  of  Z  by  cannibali¬ 
zation.  We  assume  that  for  any  value  of  X  cannibalization 
will  be  effected  in  such  a  way  that  the  maximum  value  of  Z 
for  that  X  will  be  achieved.  See  Hirsch,  Miesner  and  Boll 
[Ref.  5,  p.  336-342]  for  a  detailed  description  of  the 
cannibalization  operation  for  the  type  of  structure  being 
discussed  here.  Let  Z*  denote  the  state  of  the  system 
after  an  optimal  cannibalization,  and  let  Z*  denote  the 
state  of  the  system, with  respect  to  FSNj  under  a  policy 
of  optimal  cannibalization. 


Note  that  Z*  may  be  considered  to  be  a  function  of 
Wj  the  number  of  working  parts  of  type  j ,  since  Z^  is 


the  value  of  the  number  of  operational  unit  systems  after 
an  optimal  .cannibalization  when  the  only  part  failures  are 
at  those  loci  associated  with  part  type  j.  It  is  assumed 
that  Z^(Wj)  is  a  monotone  non-?decreasing  function  of  W ^ . 
Let  K (Wj , z)  denote  the  set  of  values  of  the  number  of 
working  parts  W .  such  that  the  value  of  Z*  is  at  least 
as  great  as  z„  i.e., 

K(W,,z)  =  {W.:  0,  <  Wj  <  Sa.  and  Z*(W ..)  i  z). 

Note  that 

K  (W.. , 0 )C  K  (W.. ,  1) C  •  • .  C K  (Wj  , S)C  ^  (W  .  , S+l)  ,  j=l , . . .  ,N; 

and  sirice  max  Z*  =  S, 

K  (Wj  /S+l)  =  <f>  . 

The  largest  value  of  z  such  that  K(Wj,z)  ^  <J>  is  S. 

Further  note  that  the  minimum  value  of  such  that  Z 

Z*  =  z  ivs  zb .  since  b .  is  the  minimum  number  of  part 
3  3  3 

type  j  needed  on  each  unit  system  if  that  unit  system  is 
to  be  operational,  in  symbols 

zbj  =  min  {w^:  WjSK(W..,z)}  . 

Let  Kex  be  an  arbitrary  set  where  ^  is  the  set  of  all 
system  state  vectors  X.  The  indicator  function  lv  is 
defined  for  all  K  and  X  as  follows 

IK(X)  =  1,  if'  XeK, 

=  0,  if  X^K  . 
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Now  let  us  use  the  symbol  {Wj  >.  zb^}  to  denote  the  set  of 
all  X  such  that  w^  >.  zbj,  i.e.,  {Xe$(:  Wj  (X)  -  zb^}. 

Hirsch,  Meisner  and  Boll  show  (using  somewhat  different 
notation)  that  if  Z*  =  min  Z*. ,  then 

j  3 


S  S  N 

Z*  =  -  I  M  =  E  n  L  s  i  , 

k=l  ^  {Wj— zb  j}  k=l  j=l  {wj  “  zbj} 

j=l 

and  the  authors  call  this  the  "representation  theorem"  [Ref. 
5,  p.  349],  If  S  =  1,  i.e.,  we  have  only  one  unit  system 
then  Z*  =  Uf  and  U  =  1  if  and  only  if  each  one  of  the 
indicators  I ,  .  ■>  is  equal  to  1.  The  structure  of 

‘"rV 

the  unit  system  in  this  model  is  thus  a  coherent  structure 
in  the  definition  of  that  term  given  by  Birnbaum,  Esary  and' 
Saunders  [Ref.  4,  p.  61].  Further  each  unit  system  may 
be  considered  to  be  a  k-out-of-n  structure  with  respect 
to  each  part  type  j ,  where  a  k-out-of-n  structure  is  one 
which  has  n  parts  of  a  single  type  and  is  operational  if 
and  only  if  at  least  k  out  of  the  n  parts  are  working. 
Also /  we  can  consider  each  unit  system  to  be  a  series  system 
composed  of  N  k-out-of-n  structures/  where  the  k's  are 
the  bj's  and  t.he  n's.  are  the  a^'s. 


B.  .CORRELATION  OF  INVENTORY  POLICIES  AND  SYSTEM  AVAILABILITY 
Consider  now  the  effect  of  inventory  policies  on  the 
expected  value  of  the  number  of  operational  unit  systems 
E(Z).  If  there  were  an  infinite  number  of  repair  parts  of 
each  type  available,  the  value  of  X  would  always  be 
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1  =  (1,1,... ,1) ,  since  it  is  assumed  that  each  part  is 
replaced  immediately  upon  failure  as  long  as  spares  are 
available,  and  the  state  of  the  system  Z  would  always 
be  S.  There  are  in  fact  various  constraints  on  the  ability 
of  the  support  unit  to  carry  inventories  of  parts.  The 
result  is  that  from  time  to  time  demands  will  occur  for 
which  no  replacement  part  is  immediately  available.  For 
any  state  vector  X,  denote  the  state  vector  after  an 
optimal  cannibalization  as  X*  and  the  value  of  Z  after 
cannibalization  as  Z*.  Then 

X*(t)  =  (X*  ^(t) X|  (t) )  ,  t  >  0, 


is  a  stochastic  process  where  for  each  fixed  t,  Xk^ 
a  random  variable  taking  the  values  zero  or  one.  Now 
each  subset  K  of  the  set  of  possible  values  of 
let  Pt(K)  denote  the  probability  that  X*  (t)  is  an 
ment  of,  the  set  K,  ie., 


is 
for 
X  , 
ele- 


Pt(k)  =  Pr{X*(t)eK},  KCX  • 


As  shown  by  Hirsch,.  Meisner  and  Boll  the  probability 
distribution  'defined  above  concentrates  all  of  its  mass  on 
the  set -of  maximum  points-  Mx,  which  are  the  possible 
vectors  X*,  the  state  vector  after  an  optimal  cannibali¬ 
zation.  Any  vector  X  in  which  the  number  of  working 
parts  of  each  type  is  given  by  the  vector  W  can  be 
transformed  into  another  vector  X’  which  has  the  same 
value  of  W.  Thus,  since  an  optimal  cannibalization  yields 
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the  maximum  value  of  the  system  structure  function  Z  for 
a  given  state  vector  X,  there  is  a  unique  value  of  the 
state  of  the  system  after  an  optimal  cannibalization  Z* 
for  each  value  of  the  vector  W  =  (W^,...,WN).  The  set  of 
all  X  corresponding  to  a  given  W  vector  is  called  an 
equivalence  class,  and  the  set  is  a  subset  of  this 

equivalence  class  where  Z  as  a  function  of  X  takes  its 
maximum  values  over  the  equivalence  class.  We  denote  each 
X  in  Mx  as  X*.  To  describe  the  variation  in  time  of 
the  system  state  function  Z,  set 

Z*(t)  =  Z  (X*  (t ) )  ,  t  0, 

and  we  note  that  the  probability  distribution  of  Z*  (t)  is 
given  by 


Pr{Z*(t)  >  z}  =  Pr{Z*  (x*  (t) )  z} 

=  Pt('X:  Z*(X)  z}  . 

Let  Wj (t)  represent  the  number  of  working  spares  of 
type  j  at  time  t.  Under  the  assumption  that  Z*  =  min  Z*, 

j  j 

we  have  immediately  from  the  representation  theorem  that 

z*!t)  =  k=i  = 

Assume  that 

1.  For  each  index  j  =  1,...,N,  the  loci  in  which 
parts  of  type  FSNj  occur  are  indistinguishable  in  their 
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effects  on  the  lifetimes  of  the  parts  installed  in  them, 
i.e.,  the  failure  rate  at  each  instant  of  a  given  part  type 
FSN j  does  not  depend  on  the  locus  in  G(FSNj)  in  which 
the  part  is  installed,  nor  on  the  particular  sequence  of 
loci  through  which  it  has  passed. 

2.  Parts  operate  independently,  i.e.,  the  lifetime  of 
a  given  part  is  not  related  to  the  lifetimes  of  any  other 
parts.  These  assumptions  make  it  reasonable  to  postulate 
that  the  joint  distribution  of  (W^ (t) , . . . ,WN (t) )  does  not 
depend  upon  the  particular  cannibalizations  involved  in  the 
process  {Z*(t),  t  >  0  and  that  the  N  stochastic  pro¬ 
cesses 

{Wx(t)  ,  t  >  0},  . . . ,  (WN(t)  ,  t  >  0} 

are  mutually  independent.  Thus  the  expected  number  of  unit 
systems  up  under  optimal  cannibalization  at  time  t, 
E,{Z*(t)),  is  given  by 

E(z*(t)>  =  jj.  Jj.  E[1<w.(t)  >  kbj)l 

S  N 

=  E  JI  Pr{W . (t)  >  kb . 
k=l  j=l  3  3 

Further, the  probability  that  Z*  (t)  equals  at  least  z 
is  given  by 

N 

Pr{Z*  (t)  z}  =  n  Pr {W .  (t)  >  zb.}. 

j=l  3  3 
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Now  let  Yj  (t)  be  a  random  variable  denoting  the  number 
of  backorders  for  part  type  FSN  j  at  time  t.  Under  the 
assumption  that  every  failed  part  is  replaced  immediately 
by  a  working  part  from  the  inventory  of  spares,  we  see  that 

Wj(t)  =  Saj  -  Yj  (t) 

and 

{Wj(t)  >  kbj}  =  {Saj  ~  Yj  (t)  kbj}  =  {Yj  (t)  <  Saj-kbj}. 

Hence 

N 

Pr{zMt)  >  k}  =  HPr{Y.(t)  ;>  Sa.  -  kb.} 

j=l  3  33 

and 

S  N 

E  (Z*  (t) )  -l  2  Pr{Y.  (t)  >  Sa.  -  kb.}. 

k=l  j=l  D  J  * 

We  assume  that  the  inventory  system  is  continuous  review 
with  a  (Qj,rj)  policy  for  each  FSN ,  where  Qj  denotes 
order  quantity  and  rj  denotes  reorder  point.  We  also 
assume  the  lead  time  demand  is  Poisson  distributed.  Hadley 
and  Whitin  [Ref.  8,  p.  184] show  that  under  a  (Q,r)  policy 
when  the  lead  time  demand  is  Poisson  with  parameter  px 
that  the  probability  that  there  are  y  backorders  at  an 
arbitrary  time  t  under  steady  conditions  may  be  derived 
as  follows.  The  inventory  position,  IP,  is  defined  as  the 
stock  on  hand  plus  stock  on  order  minus  backorders.  The 
IP  varies  between  r+1  and  r+Q,  where  r  is  defined  in 
terms  of  the  IP.  Further,  Hadley  and  Whitin  show  that  the 
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probability  that  the  IP  is  in  any  state  r+ j ,  j=l,...,Q 
is  1/Q.  Thus  if  the  inventory  system  is  in  state  r+j  at 
time  t-T,  the  probability  of  y  backorders  at  time  t 
equals  the  probability  that  y+r+j  demands  occur  during 
the  period  (t-x)  to  t  for  y>0.  The  probability  that 
y=0  is  the  probability  that  the  demand  during  (t-x,t) 
is  less  than  or  equal  to  r+j.  Hadley  and  Whitin  do  not 
make  this  distinction  for  y=0 ;  consequently,  the  formula 
they  derive  is  not  valid  for  y=0.  When  0  is  substituted 
for  y  in  the  expression  they  derive  the  resulting  value 
is  the  probability  that  y=0  and  on  hand  inventory  =  0. 

Terms  for  the  probability  that  y=0  and  on  hand  inventory 
is  greater  than  zero  are  left  out  of  their  expression.  The 
expression  Hadley  and  Whitin  derive  is 

Q  y+r+Q 

ip  (1/Q)  E  p  (y+r+j  ;px)  =  (1/Q)  2  p(u;px) 

-  '  j=l  u=y+r+l 

=  (1/Q)  [P  (y+r+l;px)  -  P  (y+r+Q+1; px) ] , 

where  p(u;px)  is  the  Poisson  probability  mass  function  with 

parameter  px  at  the  point  u  for  u  =  0,1,2,...;  and 
00 

P(u;px)  =  Z  p(v;px)  is  the  Poisson  complementary  cummulative 
v=u 

distribution  function.  The  above  expression  is  valid  y 
greater  than  zero.  For  y  =  0  , 

Q  r+j 

'  ^(0)  =  (1/Q)  E  Z  p (u; px) 

j=l  u=0 
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Q 

=  (1/Q)  2  [1  -  P(r+j+l;PT)] 

j=l 


Q 

=  1  -  (1/Q)  2  P(r+j+l;pT). 

j=l 


We  assume  that  the  pooled  demand  rate  pj  for  all 
parts  of  type  j  is  given  by 


p.  =  L  Sa.p.  , 
J  a  j  J 


where  L_  is  the  activity  level  multiplier  discussed  iii 

Section  I S  is  the  number  of  unit  systems,  a^  is  the 

number  of  applications  of  part  type  j  on  each  unit  system, 

and  p .  is  the  failure  rate  for  part  type  j  when  L 
J  ® 

is  one. 


Let  t j  be  the  procurement  lead  time  for  part  type  j. 
Demand  is  assumed  to  be  Poisson  distributed  with  parameter 


PjTj.  Thus  in  terms  of  ip ^  (y ^ ) , 


Sa.;-kb^ 

Pr{Z*(t)  2.  k}  =  ZJ  3  .  (y  . )  . 

J  y  ,  =  0  33 


Also  since  Z*(t)  =  min  Z*(t)  and  the  y.  are  assumed 

j  3  3 

to  be  independent,  the  probability  that  Z*(t)  is  equal 
to  or  greater  than  k  is  'he  probability  that  each  Z*  (t) 
is  equal  to  or  greater  than  k,  i.e., 


N  Sa^-kb. 

Pr  fc*(t)  >  k}  =  n  E-  3  i|>.  (y .)  . 

j=l  y ^=0  3  3 
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Now  the  expected  value  of  a  non-negative  integer  random 
variable  is  the  sum  from  1  to  *  of  the  complementary 
cumulative  probability  function;  thus 

CO  S 

E [Z* (t) ]  =  Z  {Pr  Z*  (t)  >  k}  =  Z  Pr{Z*  (t)  >  k}, 
k~l  k=l 

since  the  probability  that  Z*  (t;  is  greater  than  S  is 
zero.  Therefore,  the  expected  value  of  the  number  of 
operational  systems  in  terms  of;  the  probability  mass 
functions  of  the  backorders  for  each  part  type  is: 

S  N  Sa-i-kb.! 

E[Z*(t).]  =  2  n  SJ  J  il>.( y.)f 

k=l  j=l  y  j=G  3 

where: 

xp j  (y  j )  =  (l/Q)  [P  (yj+r^+l)  -  P  (yj+r^+Qj+l)  ] ,  for  y^  >0 , 

and 

Q 

(Y  • )  =  1  -  (l/Q)  2  P(r.+m+l),  for  y.  =  0. 

J  D  m=l  J  3 

The  problem  of  optimum  inventory  policies  may  be  formulated 
as  the  minimization  of  the  expected  cost  of  on  hand  inven¬ 
tory  subject  to  one  of  the  following  constraints':  first, 
that  the  expected  number  of  operating  systems  E(Z)  must 
be  equal  to  or  greater  than  some  required  fraction  of  the 
total  number, or,  second,  that  the  probability  that  at  least 
k  of  the  unit  systems  are  operational  must  be  equal  to  or 
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greater  than  a  minimum  assurance  level.  The  expressions 
used  to  relate  inventory  policies  and  these  constraints 
are  those  developed  above  for  E(Z)  and  Pr{Z  k.  k}.  Note 
that  since  the  expressions  refer  to  "steady  state"  con¬ 
ditions,  and  since  optimal  cannibalization  is  always 
performed,  the  symbol  Z  will  be  used  henceforth  in  lieu 
of  Z*  (t)  and  Z  will  indicate  the  "steady  state"  random 
variable  Z* (t) . 

It  turns  out  that  for  part  types  with  low  demand  rates 
the  optimal  policy  under  the  formulation  of  the  problem 
indicated  above  is  not  to  stock  these  low  demand  items 
at  the  direct  support  level,  but  to  order  them  from  the 
source  of  supply  as  demands  occur.  If  this  rs  the  case, 
r  will  be  set  to  -1  and  Q  to  1.  This  means  that  an 
order  is  placed  as  soon  as  a  failure  occurs.  In  this 
situation  the  distribution  of  Yj  reduces  to  the  Poisson 
distribution  with  parameter  PjTj'  since  the  inventory 
position  is  always  equal  to  zero,  and  there  is,  therefore, 
only  one  state  for  the  inventory  position.  Hence  the 
probability  that  there  are  exactly  yj  backorders  at  any 
arbitrary  time  is  p(yj;PjTj). 

It  is  interesting  to  consider  what  would  happen  if 
PjTj  were  small  for  all  N  part  types,  and  no  stock  were 
carried  at  the  direct  support  level.  Assume  that  each  part 
type  were  always  available  from  the  source  of  supply  with 
a  procurement  lead  time  of  ij,  a  constant,  which  might  be 
different  for  different  part  types.  Now,  if  no  parts  are 
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carried  at  direct  support  level,  but  are  ordered  only  when 
-demands  occur,  i.e.,  r^  =  -1,  Q-  “  1  for  all  j,  then 

<Pr{Z  >.  k}  =  II  E^  ^  p(yi;p.t.)/ 
j=l  y ^=0  3  3  3 

and 

S  N  Sa-i-kb . 

E  (Z)  •-=  E  II  EJ  3  p(y.;p.T.). 
k=l  j=l  y j=0  3  3  3 

The  values  of  the  above  functions,  can  be  used  to  estimate 
the  lower  bounds  for  E(Z)  and  Pr{Z  k},  when  demand 
rates  and  procurement  lead  times  are  p\.,T.,  j=l,...,N. 

J 
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III.  OPTIMIZATION  TECHNIQUE  AND  NUMERICAL  ANALYSIS 


A.  THE  COST  FUNCTION 

The  objective  function-  to  be  minimized  is  a  cost 
function  of  the  following  form: 


E (C)  =  Z  C.D(Q.,r.)  , 
j=l  3  33 


where 


E (C)  =  Expected  cost  of  on  hand  inventory  in 

dollars/  cubic  feet/  pounds  or  some 
other  measure  of  a  resource  in  short 
supply, 

Cj  =  The  unit  cost  of  part  type  j, 

D(Q. ,r.)  -  The  expected  on  hand  inventory  of  part 
■f  3  type  j ,  given  order  quantity  ,Q  . 

and  reorder  point  r^.  3 


If  E (C)  is  in  dollars,  the  above  function  gives  the 
expected  amount  of  funds  tied  up  in  repair  parts  inventory 
for  the  system  at  the  direct  support  echelon.  If  E(C) 
were  in  pounds  the  function  might  represent  the  expected 
amount  of  load  carrying  capacity  of  the  maintenance  used  up 
by  the  repair  parts  for  the  system  in  question. 

The  value  of  D(Qj,rj)  is  given  by 

D(Q.,r.)  =■  (Qj+l)/2  +  r .  -  p.T.  +  B(Q..,r.), 

where 


B(Q,,r.)  =  Z  y.iMy,;Q.,r.) 
J  J  „  _Q  J  J  J  J  J 

X  j 
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is  the  expected  number  of  backorders,  given  the  inventory 
policy  (Qj,rj).  Hadley  and  Whitin  [Ref.  8,  p.  184]  show 
that 

00 

B (Q,r)  =  (1/Q)  E  y  [P<(y+r+l; px)  -  P (y+r+Q+1; pt)  ] . 
y=o 

B.  OPTIMIZATION  TECHNIQUE 

1.  Minimizing  Cost  Subject  to  a  Required  Expected 
Number  of  Unit  Systems  Up 

Suppose  it  is  required  that  the  expected  number  of 
unit  systems  up  must  be  equal  to  or  greater  than  A_S, 

o 

where  A  is  the  required  availability  with  respect  to 
supply  support,  0  £  A  £.1,  and  S  is  the  number  of  unit 

S 

systems.  Further,  suppose  it  is  desired  to  minimize  the 
cost  of  expected  on  hand  inventory  needed  to  achieve 
E  (Z)  >.  A  S.  The  problem  may  then  be  stated  as 

b 

N 

Minir  ize  E  C.D(Q.,r.) 

(r1;...,rN)  j=l  3  3  3 

(Q]_ ,  •  •  •  ) 

subject  to  E  (Z)  >_  AgS. 

Finding  the  optimal  solution  to  this  problem  for 
(Q j ,r> j )*',  j  =  1,...,N  is  a  formidable  task  due  to  the  com¬ 
plexity  of  the  cost  function  and  the  function  for  the 
expected  value-  of  Z  in  terms  of  the  Q .  and  r^,-  To 

J  J 

simplify  the  problem  we  have  chosen  to  set  the  values  of 
Qj  by  an  operating  level  of  supply  in  terms  of  days  or 
months  of  supply.  Thus  we  set 
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Qj  =  [pjL0  +  .5],  if  [PjL0  +  .5]  >0, 

=  1,  otherwise, 

where:  [u]  =  The  greatest  integer  equal  to  or  less  than  u, 

L  =  The  operating  level  of  supply,  e.g.,  one  month 
of  supply. 

The  problem  now  becomes:  given  the  values  of  Q..  determined 
as  indicated  above,  find  the  vector  (r^,...,rN)  which 
minimizes  the  cost  function  subject  to  the  required  level 
of  availability.  We  shall  restrict  the  r^  to  be  >  -  1. 
After  computing  the  values,  the  next  step  in  the  pro¬ 

cedure  is  to  compute  a  table  of  values  of  i|k  (y ^ ;  r  j  =  -1) 
for  all  j .  This  table  is  useful  in  that  the  values  of 
may  be  found  for  other  values  of  r^  as 

follows :. 

r  j+1 

(0;  r.)  =  E  iMu;  -1),  for  y.  =  0, 

J  J  u=0  J  J 

rj}  =  *j  (yj+rj+1;  ”1)'  for  yj>0* 

The  number  of  values  needed  for  each  depend  upon  the 

demand  rate  and  the  accuracy  desired  in  the  computations. 
After  the  table  of  values  is  computed  the  next  step  is 

4 

to  find  a  set  (rj, . .  ,r*)  of  r^  such  that  E(Zj),  the 
expected  value  of  Z  with  respect  to  part  type  j,  is 
equal  to  or  greater  than  A  S,  and  such  that  Pr(Z.=S;r.) 

s  J  J 

is  concave  for  r^  —  rj*  Each  r^  must  be  large  enough 
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to  insure  E(,Zj)  >.  AgS  since  E(Zj)  gives  the  expected 

number  of  operating  systems  when  every  other  type  except  j 

is  always  in  supply.  E(Z)  can  not  be  greater  than  E(Z^) 

because  it  is  assumed  that  at  any  arbitrary  time  Z  -  m^n  Z . . 

3  3 

Now  E  (Z j )  is  given  by 


E  (Z  j ) 


S  Sa --kb  • 


k=l  y.=0 


r.) 


Using  this  expression  find  the  minimum  r..  such  that 

E (Z  . )  2.  A  S.  (See  Appendix  A  for  a  discussion  of  the  need 
J  S 

to  require  each  Pr(Z.  >.  k;  r.)  to  be  concave  for  r.>r^.) 

3  3  3  3 

After  the  set  of  initial  values  of  the  r..  is  found,  a 
marginal  analysis  technique  similar  to  the  approach  of  Karr 
and  Geisltr  [Ref.  9]  and  to  one  form  of  optimal  redundancy 
algorithm  of  Barlow  and  Prochan  [Ref.  3,  p.  166]  is  used  to 
find  at  each  iteration  the  part  type  which  yields  the 
greatest  increase  in  E(Z)  for  the  increase  in  the  cost  of 
the  expected  on  hand  inventory  for  that  part  type.  On  the 
first  iteration  the  ratio 


__  E(Z;rj+l)  -  E  (Z?r^) 
Aj  -  CjlDj1rj+l)  -  Djl2j: 


is  computed  for  each  part  type,  where 


S  Pr (Z . >k;r .+1) 

E(Z;r.+l)  =  S  Pr(Z^k)  =  j-r- 


Pr(Z.>k;  r.)  * 
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The  maximum  over  j  of  A^  is  determined  and  the  r^ 
corresponding  to  that  maximum  is  increased  by  one.  E(Z) 
is  set  equal  to  E{Z;rj+l)  corresponding  to  this  j.  If 
the  new  value  of  E(Z)  is  equal  to  or  greater  than  AS, 
the  procedure  stops;  if  not,  a  new  set  of  A^  is  computed 
and  the  r^  corresponding  to  the  maximum  over  j  of  the 
new  Aj  is  increased  by  one.  The  procedure  continues  in 
this  fashion  until  the  value  of  E(Z)  is  equal  to  or 
greater  than  AsS. 

2 .  Minimizing  Cost  Subject  to  a  Required  Probability 
that  at  Least  k  Unit  Systems  are  Up 

Suppose  that  it  is  required  that  P{z  >.  k}  >. 

P  .^(k),  and  that  it  is  desired  to  minimize  the  cost  of 
min 

the  expected  on  hand  inventory  needed  to  achieve  this 
probability.  The  general  procedure  for  finding  the  optimum 
vector  (r^,...,rN)  is  the  same  as  the  E(Z)  :>.  AsS  form 
of  the  constraint  with  the  following  expections: 

a.  The  initial  r ^  values  (rj,...,r*)  are  the 

minimum  values  such  that  Pr(Z.  2.  k)  >.  P  .  (k)  and 

3  nun 

Pr(Z.  >.  k.r.)  is  concave  in  r.  for  r.  >.  r*. 

3  ]  ]  D  3  3 

b.  At  each  iteration  Aj  is  computed  as 

Pr{Z  >.  k;r ^+1}  -  Pr{Z  >.  k;  r^} 

Aj  =  C.[D(r.+l)  -  D(rj)] 

As  in  the  first  form  of  the  constraint  the  rj  associated 
with  the  maximum  A^  is  increased  by  one  at  each  iteration 
until  Pr(Z  >.  k)  >.  P  .  (k) .  Since  Pr{Z  k.  k}  is  just  one 
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term  of  the  sum  over  k  in  the  expression-  for  E(Z)  the 
computational  effort  at  each  iteration  is  much  reduced  for 
this  form  of  the  constraint,  particularaly  if  the  number 
of  unit  systems  S  is  large.. 

C.  NUMERICAL  ANALYSIS  AND  DEVELOPMENT  OF  FORTRAN  PROGRAMS 
Two  FORTRAN  programs,  one  for  each  form  of  the  constraint 
have  been  developed  and  run  for  several  sample  problems  on 
the  Naval  Postgraduate  School's  IBM  360  computer.  The  pro¬ 
grams  are  very  similar,  and  both  use  identical  versions  of 
the  following  subprograms:  PSITAB,  which  computes  a  table 
of  values  of  ip^(y;r^.  =  -1)  for  j  =  1,...,N;  PRZJK,  a 
function  subprogram  which  computes  Pr{Z.  >.k;r.};  EBO,  a 

•J  J 

function  which  computes  the  expected  number  of  backorders 
for  part  type  j  give  that  the  reorder  point  is  r ^ ; 

PPT,  a  function  which  computes  individual  terms  of  the 
Poisson  probability  mass  function;  and  PCUMT,  which  computes 
complementary  cumulative  terms  of  the  Poisson  distribution. 

With  the  IBM  360  it .was  found  necessary  to  do  the 
arithmetic-  in  double  precision  because  the  single  precision 
round-off  errors  for  repeated  multiplications  were  excessive. 
The  table  of  (y ; r j  =  -1)  values  was,  however,  stored 

as  a  single  precision  array  to  conserve  storage  space.  When 
values  from  this  table  were  needed  they  were  converted  to 
double  precision  with  the  standard  function  DBLE . 

Subroutine  PSITAB  stores  the  computed  values  of 
(y 7 r j  =  "D  in  an  array  called  PI(J,K),  where 
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PI  { J /  K)  =  'i'j  (k-l;rj  =  -1) . 


The  array  size  is  (N,100) ,  where  N  is  the  number  of 
different  part  types;  however,  for  each  J  only  values 
of  PI ( J ,K)  for  K  less  than  or  equal  to  KMAX(J)  are 
computed,  where 


KMAX(J)  =  min  JlOO,  largest  K  such  thati 

PI ( J ,K)  <  10“12 


PI,(J,K)‘  for  K  greater  than  KMAX(JT)  are  set  to  zero. 

For  part  types  with  low  demand  rates  the  number  of  values 
of  PI ( J,K)  computed  is  much  less  than  100.  The  array 
size  used  in  the  programs  permits  handling  demand  rates 
up  to  about  50  without  significant  truncation  of  the  distri¬ 
bution  of  Yj . 

1.  Program  EZMIN 

Program  EZMIN  computes  the  optimum  r  vector  sub¬ 
ject  to  a  required  minimum  expected  value  of  the  number 
unit  systems  up.  The  main  program  reads  the  input  data, 
writes  the  system  parameters  and  the  required  availability, 
and  calls  subroutines  PSITAB,  INITAL,  and  OPTIMR.  Sub¬ 
routine  PSITAB  computes  the  PI(J,K)  array  as  discussed 
above.  Subroutines  INITAL  computes  an  initial  set  of  r^ 
values,  r£,...,r*,  each  of  which  satisfies  the  following 
conditions:  Pr(Zj  =  Syr^)  a  concave  function  of  r^  for 

r.  >  r*; 

3  ~  3 ' 


E  (Z  . )  =  E  Pr (Z  .  >  k?r*)  >  AS. 
J  \r —  'I  J  J  S 
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Subroutine  OPTIMR  uses  the  marginal  analysis  technique 
previously  discussed  to  find  the  optimal  r  vector. 

2.  Program  PKSMIN 

Program  PKSMIN  computes  the  optimum  r  vector 
subject  to  a  required  assurance  that  at  least  k  out  of 
S  unit  systems  are  operational.  This  program  differs  from 
program  EZMIN  only  in  the  following  respects.  PKSMIN  Main 
calls  subroutine  INTLZ  instead  of  INITAL  and  subroutine 
OPTMZ  instead  of  OPTIMR.  Subroutine  INTLZ  finds  an 
initial  set  of  r^  values  r|,...,r^  such<  that  Pr(Zj>k;rj) 
is  concave  for  r.  >  r*  and  such  that  Pr(Z.  >  k;r*)  >.  P 
Subroutine  OPTMZ  computes  the  optimum  r  vector  subject  to 
the  required  assurance  that  Z  be  equal  to  or  greater  than 
k. 

3 .  Flow  Charts 

Flow  charts  for  program  EZMIN  are  presented  in 
Appendix  B.  Flow  charts  for  program  PKSMIN  are  presented 
in  Appendix  C. 
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IV.  EXAMPLE  PROBLEMS  AND  DISCUSSION  OF  RESULTS 


A.  DESCRIPTION  OF  THE  EXAMPLE  SYSTEM 

The  example  system  consists  of  50  10S-mm  recoilles 
rifles,  model  M40A1.  Several  factors  made  this  system  a 
convenient  example.  First,  the  total  number  of  repair 
parts  is  small  enough  so  that  the  FORTRAN  programs  developed 
for  solution  of  the  two  alternate  formulations  of  the  prob¬ 
lem  can  be  used  without  modification  on  the  IBM  360  G-level 
computer  available  at  the  Naval  Postgraduate  School. 

Second,  Department  of  the  Army  Technical  Manual  9-1015-221-35 
[Ref.  11]  contains  data  from  which  failure  rates  may  be 
estimated.  Third,  actual  demand  data  for  a  six  month 
period  for  62  weapons  used  in  training  heavy  weapons  infan¬ 
trymen  at  Fort  Ord,  California,  were  made  available  to  'the 
author  by  the  Fort  Ord  Post  Maintenance  Section.  Finally, 
price  data  were  available  from  the  microfilm  Army  Master 
Data  File  Selected  Management  Data  File  [Ref.  11], 

The  parts  list  in  TM  9-1050-221-35  contains  287  different 
Federal  Stock  Numbers  (FSN) .  Of  these  only  the  I59  FSN 
which  were  indicated  as.  being  combat  essei._ial  were  used 
for  the  sample  problems.  By  limiting  consideration  to 
only  combat  essential  parts  we  insure  that  some  of  the 
assumptions  of  the  model  are  more  nearly  satisfied  than  if 
all  part  types  were  considered.  For  example >  the  model 
assumes  that  for  each  part  type  j  there  is  a  positive 
number  bj  which  is  the  minimum  number  of  working  parts  of 
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type  j  needed  for  the  unit  system  to  be  capable  of  per¬ 
forming  its  function,  A  part  would  not  be  combat  essential 
if  bj  were  zero.  Thus  the  assumption  that  the  unit  system 
is  a  series  of  k-out-of-n  structures  and  that  the  possible 
states  are  zero  or  one  is  more  nearly  true  if  the  unit 
system  consists  only  of  combat  essential  parts. 

The  Federal  Stock  Number,  the  number  of  applications 
per  unit  system,  and  data  for  estimating  failure  rates 
were  obtained  from  columns  2 ,  4 ,  andi  6 ,  respectively  of 
Section  2,  Appendix  B,  Ref.  10. 

B.  SAMPLE  PROBLEM  1  -  REQUIRED  EXPECTED  NUMBER  OF  UNIT 

SYSTEMS  UP  CONSTRAINT 

The  first  sample  problem  was  computed  for  the  required 
expected  number  of  -unit  systems  up  form  of  the  constraint. 
Table  I  lists  a  summary  of  the  input  data  and  results  for 
this  problem. 


TABLE  I 

SAMPLE  PROBLEM  1  -  SUMMARY  OF  INPUT  DATA  AND  COMPUTED  RESULTS 


INPUT 

DATA 

NUMBER  OF  UNIT  SYSTEMS 

50 

NUMBER  OF  PART  TYPES 

159 

REQUIRED  AVAILABILITY 

0.95 

OPERATING  LEVEL 

TTo 

ACTIVITY  LEVEL 

o 

• 

rH 

COMPUTED 

RESULTS 

EXPECTED  NUMBER  OF 
OPERATIONAL  UNIT  SYSTEMS 

47.58 

EXPECTED  COST  OF 

ON  HAND  INVENTORY 

$418.04 

NUMBER  OF  ITERATIONS  OF 
MARGINAL  ANALYSIS  PROCEDURE 

295 
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The  above  data  show  that  for  the  system  consisting  of 
50  weapons  with  159  part  types  and  a  required-  availability 
of  0.95  program  E2MIN  computed  an  r  vector  with  an 
associated  expected  cost  of  on  hand  inventory  of  $418.04. 

The  computed  value  of  E(Z)  is  47.58,  which  indicates  that 
the  marginal  analysis  technique  overshot  the  required 
expected  number  of  systems  up  by  .08,  since  the  required 
expected  number  was  50  x  0.95  =  47.50. 

Table  II  lists  the  input  parameters  and  computed 
results  for  each  FSN.  Parts  are  listed  in  Federal  Item 
Identification  Number  (FUN)  sequence,  because  they  are 
listed  this  way  on  the  Army  Master  Data  File  [Ref.  11] . 

The  FUN  is  the  last  seven  digits  of  the  Federal  Stock 
Number.  The  data  cards  were  sorted  into  this,  sequence  to 
facilitate  determination  of  unit  prices.  The  columns  of 
Table  II  contain  the  following:  j,  the  sequence  number; 
FSN,  the  Federal  Stock  Number;  A  =  a ^ ,  the  number  of 
applications  of  part  type  j  on  each  unit  system;  B  =  b^, 
the  minimum  number  of  part  type  j  needed  for  the  unit 
system  to  be  up;  C  =  C.,  the  unit  cost  of  part  type  j; 
RHOHAT  =  p . ,  the  failure  rate  for  part  type  j ;  RHO  =  p .  = 

J  J 

Sa^Pj,  the  pooled  demand  rate  for  all  parts  of  type  j  in 
the  system;  TAU  =  tj,  the  procurement  lead  time  for  part 
type  j;  Q  =  Q j ,  the  order  quantity;  R  =  r ^ ,  the  reorder 
point;  EB  =  B(Qj,rj),  the  expected  backorders  for  part  type 
j;  and  EOH  =  D(Qj,rj),  the  expected  on  hand  inventory  for 
part  type  j . 


TABLE  II 

SAMPLE  PROBLEM  1  rr  REQUIRED  EXPECTED  NUMBER  UNIT  SYSTEMS  UP  CONSTRAINT  - 
INPUT  PARMETERS  AMD  COMPUTED  VALUES  ,OF  ORDER  QUANTITIES,  REORDER  POINTS 
EXPECTED  BACKORDERS  AMD  EXPECTED  ON  HAND  QUANTITIES  FOR  ALL  PART  TYPES 


ID 

fA 

ID 

CM  CM 

CM 

D-  CM 

jt 

CD 
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C.  SAMPLE  PROBLEM  2  -  REQUIRED  ASSURANCE  AT  LEAST  K  UNIT 
SYSTEMS  UP  CONSTRAINT 

The  second  sample  problem  was  computed  for  the  required 
assurance  at  least  k  unit  systems  are  up  form  of  the 
constraint.  Table  III  lists  a  summary  of  the  input  data 
and  results  for  this  problem. 


TABLE  III 

SAMPLE  PROBLEM  2  -  SUMMARY  OF  INPUT  DATA  AND  COMPUTED  RESULTS 


INPUT- 

DATA 

NUMBER  OF  UNIT  SYSTEMS 

50 

REQUIRED  NUMBER  OF  UNIT 
SYSTEMS  UP 

47 

REQUIRED  ASSURANCE 

0.90 

OPERATING  LEVEL 

1.0 

ACTIVITY  LEVEL 

1.0 

COMPUTED 

RESULTS 

PROBABILITY  AT  LEAST  K 

UNIT  SYSTEMS  ARE  UP 

0.91 

EXPECTED  NUMBER  OF 

UNIT  SYSTEMS  UP 

47.46 

EXPECTED  COST  OF  ON 

HAND  INVENTORY 

$387.88 

NUMBER  OF  ITERATIONS  OF 
MARGINAL  ANALYSIS  PROCEDURE 

348 

The  above  data  show  that  for  the  same  system  considered 
in  Problem  1  with  a  required  assurance  of  0.90  that  at  least 
47  unit  systems  are  up  program  PKSMIN  computed  an  r  vector 
with  an  associated  expected  cost  of  on  hand  inventory  of 
$387.88.  The  program  overshot  the  required  probability 
by  0.01  since  the  computed  probability  is  0.91  and  the  re¬ 
quirement  was  0.90.  After  the  r  vector  which  met  the 
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constraint  was  computed  the  program  computed  the  resulting 
value  of  E(Z),  the  expected  number  of  unit  systems  up. 

It  is  interesting  to  note  that  the  expected  number  of  unit 
systems  up  very  nearly  meets  the  requirement  of  Problem  1, 
and  the  expected  cost  of  on  hand  inventory  is  also  about 
the  same  as  for  Problem  1.  The  computer  execution  time  for 
Problem  2  was  about  19  seconds  compared  to  slightly  over 
four  minutes  for  -Problem  1. 

Table  IV  lists  the  input  paramaters  and  computed  results 
for  each  FSN.  The  column  headings  are  the  same  as  for 
Froblem  1. 

D.  EXACTED  COST  OP  ON  HAND  INVENTORY  AS-  A  FUNCTION  OF 

REQUIRED  AVAILABILITY 

1.  Expected  Cost  Versus  Required  Expected  Number  of 
Unit  Systems  Up 

Program  EZMIN  was  run  for  ten  values  of  A  =  0.90, 

s 

.  ..,0.99  in  order  to  develop  a  functional  relationship 
between  the  first  form  of  the  availability  requirement  and 
the  expected  cost  of  on  hand  inventory.  Table  V  summaries 
the  results.  In  addition  to  the  expected  cost  of  on  hand 
inventory  the  table  lists  the  number  of  iterations  of  the 
marginal  analysis  procedure  required  to  find  a  solution  and 
the  number  of  reorder  points  which  were  less  than-  the  mean 
lead  time  demand  for  each  level  of  Ag.  A  reorder  point 
less  than  the  mean  lead  time  demand  is  equivalent  to  a 
negative  safety  level.  For  many  of  these  the  computed  order 
quantity  is  one  and  the  computed  reorder  point  is  minus  one. 


SAMPLE  PROBLEM  2  -  REQUIRED  ASSURANCE  AT  LEAST  K  UNIT  SYSTEMS  UP  CONSTRAINT 
INPUT  PARAMETERS  AND  COMPUTED  VALUES  OF  ORDER  QUANTITIES,  REORDER  POINTS, 
EXPECTED  BACKORDERS  AND  EXPECTED  ON  HAND  QUANTITIES  FOR  ALL  PART  TYPES 


i 


UD 

fA 

CM 

CM 

CM 

in 

ID 

CO 

fA 

CO 

00 

CO 

cn 

00 

-3“ 

CO 

CM 

00 

o 

o 

00 

fH 

CO 

CO 

o 

CO 

o 

CO 

co 

o 

in 

X 

o 

00 

r-i 

o 

o 

co 

O 

o 

o 

4 

4 

cn 

rH 

.3 

CO 

o 

tA 

o 

o 

3 

o 

co 

CD 

cn 

cn 

CM 

cn 

m 

O 

o 

o 

UD 

CM 

fA 

CM 

fA 

o 

o 

o 

CM 

o 

c 

o 

o 

CM 

00 

00 

4 

fA 

in 

CM 

O 

o 

o 

O 

o  o 

rH 

o 

-3- 

rH 

o 

c 

o 

CO 

m 

o 

c 

o 

CM 

CM 

CM 

in 

ID 

jt 

co 

o 

o 

o 

tA 

o 

ID 

00 

CD 

CT 

oo 

o 

4 

o 

co 

CM 

o 

o 

o 

00 

o 

o 

00 

rH 

ID 

ID 

o 

o 

o 

c 

o 

ID 

o 

CO 

ID 

o 

o 

in 

o 

o 

o 

00 

3 

co 

in 

rH 

O 

o 

ID 

O 

o 

o 

o 

o 

o 

o 

4 

-d 

cn 

n 

Jt 

o 

ID 

o 

CQ 

o 

o 

- t 

f^* 

a- 

o 

o 

o 

O 

O 

O 

o 

CO 

00 

•3- 

o 

CM 

CD 

o 

rH 

O 

o 

00 

in 

in 

o 

m 

o 

o 

o 

O 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

O 

O 

o 

O 

o 

o 

o 

CM 

rH 

CO 

CN 

rH 

rH 

rH 

HK>0 

ONtA 

CM  fH 

1 

1 

0 

rH 

0 

l 

0 

5 

1 

1 

i 

rH 

1 

1 

1 

1  rH 

tH 

1 

i  i 

1 

1  1 

1 

c 

rH 

4 

rH 

rH 

4 

I 

I 

4 

1 

1 

I 

I 

I 

rH 

1 

1 

1 

fA  H  H  fA 

rH 

o 

o  o  o 

o  o  o 

o  o  o 

o  o  o 

o  o  o 

o  o  o 

c 

o  o  o 

o  o 

o 

o  o  o 

coo 

o  o  o 

o  o  c 

o  o  o 

o  o  o  o 

o  o  o 

C2  o 

H 

rH 

1. 

1. 

1. 

H  H  H 

1. 

1. 

1. 

1. 

1. 

1. 

1. 

1. 

1. 

1. 

i. 

1. 

rH 

»H  rH 

rH 

rH  rH 

C 

o  a  o 

o  o  o 

o  o  o 

coo 

o  o  o 

o  o  o  o 

oo  o 

o  o 

O 

o 

o  o  o 

in  o  o 

o  o  o 

o  o  o  o  o  o 

o  o  o 

O 

o  o  o 

o  c 

a 

in  jn 

•d  -d  CO 

ID  ID  00 

id  r*  oo 

Mdd 

CMdSN 

oo  r**  oo 

cn  in 

o: 

• 

•  •  • 

•  •  • 

•  •  • 

•  •  • 

• 

• 

•  • 

• 

•  • 

o 

0 

0 

n 

0 

1 

3 

3 

5 

0 

3 

0 

0 

0 

0 

0 

0 

0 

I 

0 

CM  O 

o 

fA  rH 

H 

b 

o  o  o 

o  o  o 

o  o  o 

o  o  o 

o  o  o 

o  o  o 

O 

o  o 

o 

o  o 

< 

o 

o  o  o 

o  a  o 

o  o  o 

o  o  o 

o  o  o 

o  o  o 

o 

o  <0 

o 

o  o 

X 

4 

cn  oo 

IA  00  CM 

CM  CO  CO 

CM  Jd  CM 

CM  00  00 

4  4  4  4 

CO  -d 

co 

m  in 

o 

rH 

ft  o  o 

O  CM  S 

Mn  H 

NHO 

o  o  o 

o  r-C  I— 1 

rH 

in  rH 

rH 

o  o 

X 

o 

o  o  o 

o  o  o 

o  o  o 

o  o  o 

o  o  o 

O  o  o 

-O 

o  o 

o 

o  o 

a: 

• 

•  •  • 

•  ♦  • 

•  •  • 

•  •  • 

•  •  • 

•  •  • 

• 

•  • 

• 

•  • 

C30000  0  00000  oo  oooooooooooo 


UUnONNlAHHfflHlOtNNlOHOClNHHJlDKlMlD 
COf'.d.dVOOON'IO  NO«lDCO!flW4  O'Jt  LTV  oo  c  o  o 


OOlOOCNOOOOOOOOOJ- 


o  o 


'jHOWOO 

f- 


CO  HiOHNfOHHHOI 


00  00  H  r-t  i-C  eg  >'1 


CO 


< 

rH 

ID 

rH 

CM 

fA 

rH 

rH 

cm 

rH 

rH 

rH 

CO 

00 

ri 

rH 

rH 

CM 

rH 

rH 

rH 

rH 

rH 

4 

rH 

UD 

in 

00 

rH 

rH 

•d 

r>. 

in 

cn 

m 

4 

fA 

in 

ID 

00 

o 

rH 

in 

O 

r** 

ID 

O* 

00 

00 

UD 

rn 

rH 

ID 

•4* 

-d 

cn 

4“ 

ID 

h* 

cn 

cn 

in 

o 

o 

.4 

ID 

CM 

fA 

4 

ID 

n* 

h* 

CO 

rH 

rH 

CO 

OO 

CM 

CM 

ID 

o 

o 

o 

o 

o 

CM 

in 

in 

fA 

ID 

o 

in 

rH 

rH 

rH 

rH 

O 

CO 

C» 

O' 

fA 

in 

in 

ID 

ID 

ID 

ID 

CD 

CD 

cn 

o 

o 

in 

00 

rv 

cn 

UD 

DO 

00 

00 

00 

cn 

CM 

z 

fA 

fA 

f- 

cn 

o 

CM 

00 

00 

00 

00 

00 

cn 

00 

00 

o 

4* 

UD 

rH 

in 

cn 

cn 

cn 

cn 

in 

cn 

CO 

rH 

rH 

4 

4* 

3- 

in 

in 

in 

in 

in 

in 

in 

ID 

ID 

o 

4 

in 

1^ 

00 

rH 

rH 

rH 

rH 

CM 

ID 

u. 

o 

O 

O 

O 

O 

o 

o 

o 

o 

O' 

o 

o 

o 

o 

rH 

fl 

rH 

rH 

rH 

CM 

CM 

CM 

CM 

CM 

CM 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

m 

in 

O 

O 

O 

O 

O 

in 

in 

in 

in 

o 

in 

o 

o 

o 

tH 

rH 

rH 

fH 

rH 

rH 

rH 

rH 

rH 

o 

o 

rH 

rH 

rH 

fA 

rH 

rH 

rH 

rH 

rH 

o 

o 

fA 

m 

o 

o 

O 

m 

fA 

KA 

in 

fA 

m 

fA 

m 

m 

rH 

tH 

rH 

K. 

rH 

O 

O 

o 

fA 

H. 

in 

in 

rH 

fH 

fH 

in 

in 

in 

in 

in 

in 

in 

in 

in 

fA 

fA 

m 

in 

fA 

rH 

rH 

rH 

pH 

IA 

in 

-3 

1 

2 

3 

4 

5 

6 

7 

8 
9 
0 

rH  CM 

fA 

4 

in 

lO  N 

00 

cn 

O 

•H 

CM 

fA 

4 

in 

*H 

rH  rH 

rH 

rH 

rH 

H  H 

rH 

*H 

04 

cm 

CM 

CM 

CM 

CM 

X 

o 

UJ 


at  CM  CO  O  J  CM  00  lO  ©  ©  ©  ©  tO  tO  ©  OOCOONOO  Ot  MN 

00  00  CM  00  00  HOUjmCONW  OHO  O  00  O  rH  CM  ©  f~~  ©  CM  © 

©  rH  ©  ©  ©  O  j-cf  JK\IOO  OBS  JfCfOOOOlN  OOO 

MlOONHOOOM^NctKlKHOOOJlOOOWI^NNNffOCOHIft 
*•••  *•#••«•••••«•«••••••••••** 

ONOlA^OOOHOOOOI^OOOOOONHO^HtfOtAON 


Ct(MOMOJONCaiO|ClK>Oi010IOtOMOO»OOONOOOOiaPN 
OOOOONCOMOHOlOOlOOMflOOHDOOOWOHNCnNOWNOI 
O  WHOWIflK\OOJ4Jf(nt0OO01111SOO4ieCOJtNOOOO 
05  NOJfOONJtmOlflNWWOCTlOHOCtrfOONOOIrtNOOOO 

UJ  . . . . 

OOOOOOOOOOOOOOOCMOOOOOOOOOOOOCMO 


CO  HKlHNcfHHOHOOOOJtHHOOHHWHHCfHcfHctHO 
I  I  I  I  I  I  rH  I  I  rH 


S- 


C 1  HHrllrtKIrlrlHrlHHNHrlHICi 


N  H  rl  rl 


Lfl 


at  (O 


3 

o 

< 

• 

H 

fH. 

O 

1 

O 

o 

I 

St 

G 

Of 

'  * 

UJ 

n 

z 

•— 

H- 

o 

z 

< 

o 

o 

X 

o 

o 

rH 

X 

O 

1 

cc 

• 

o 

> 

— 

UJ 

o 

-J 

CM 

EC 

o 

• 

< 

o 

H* 

oooooooooooeoooooooooooooooooo 
ooooooooooooocooooooooooooooo 

HHHHHHHHHHHHHHHrlrlrHHHrlHHHHHHHH 

ooooooooooooooooooooooooooooo 
oooooooooocooooooooooootnooooo 

cfd'OOOIctafNlNcteooOJ'OlOf'J’d’cfONinKlNJfJMJtNO 
HO^NHOHOHOWOHCNIOOOONOHHOtAOHclW 


OOOOOOOOOOOOOOCOOOOOOOOOOO  OO O 

ooooooooooooocooooooooooooooo 

MC010«3C0C03.3C010OC00K1O30033O3UHj0WNN<i0.3O 
CMOOimNONH  CM.  rH  NNN0HH0D03HHNHN0NHN 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOrH 

OOOOCOOOOOOOOOOOOOO  O  O  OOOOOOOO 


wioa\N(aoocoo)oin(J)Ci3oiM[jHnoHfUscooNiciHwd’iciN 


OKHfWlflOONOHN3IBOOONOCO 
CM  at- 00  CM.  CM  rH  r-l 

© 


irHOOOrHOrHOOrH 


' 

j? 


00 

CM 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

CM 

rH 

ri' 

ID 

to 

rH 

rH 

CM 

CM 

rH 

rH 

CM 

rH 

rH 

< 

iM 

rH 

rH 

rH 

rH 

t-l 

rH 

rH 

rH 

rH 

rH 

CM 

rH 

rH 

ID 

to 

rH 

rH 

CM 

CM 

rH 

rH 

CM 

rH 

rH 

Cl 

in 

«o 

Cl 

O 

rH 

•st 

00 

to 

ID 

to 

r*^ 

00 

m 

m 

to 

CM 

in 

© 

O 

CM 

CM 

CM 

rH 

IM. 

Cl 

to 

OC 

00 

Cl 

Cl 

in 

rv 

rH 

CM 

CM 

CM 

CM 

rH 

m 

to 

O 

© 

at 

© 

© 

00 

© 

© 

© 

© 

IM 

TA 

co 

© 

© 

r- 

fA 

ID 

*=T 

00 

00 

00 

O 

rH 

in 

cn 

at 

CM 

CM 

O 

sn 

Cm 

© 

CM 

f- 

in 

m 

in 

in 

m 

o 

00 

rH 

fA 

© 

fA 

fA 

cn 

00 

© 

cn 

tM 

in 

© 

© 

fA 

© 

CM 

rH 

z 

to 

oa 

o 

o 

o 

O 

in 

rH 

ID 

© 

in 

in 

in 

o 

cn 

n. 

tn 

at 

in 

© 

rH 

fA 

© 

tM 

© 

CO 

Cl 

o 

o 

o 

O 

o 

CM 

© 

-if 

at 

in 

00 

o 

•=t 

fs 

o 

© 

r-l 

at 

00 

O 

It¬ 

CM 

CM 

fA 

© 

© 

© 

fA 

N> 

© 

m 

© 

fA 

fA 

fA 

m 

at 

•St 

at 

© 

© 

© 

in 

© 

© 

<£> 

O 

O 

in 

in 

in 

in 

in 

«n 

in 

in 

in 

in 

in 

o 

in 

r> 

in 

m 

© 

© 

© 

in 

© 

© 

© 

CM 

?o 

rH 

rH 

rH 

rH 

rH 

o 

© 

rH 

o 

o 

o 

fA 

o 

rH 

rH 

rH 

rH 

.rH 

O 

rH 

rH 

r-l 

© 

© 

© 

O 

O’ 

© 

O 

C* 

© 

© 

o 

o 

o 

© 

m 

fA 

tn 

O 

in 

© 

© 

O 

fA 

© 

© 

O 

in 

in 

rH 

rH 

rH 

ri 

rH 

rH 

in 

rH 

rH 

rH 

rH 

in 

m 

m 

rH 

tn 

rH 

© 

rH 

in 

© 

rH 

rH 

'D 

to 

fM 

CO 

CT1 

O 

rH 

CM 

K> 

at 

in 

to 

CO 

cn 

b 

rH 

CM 

tn 

at 

© 

© 

r*- 

oo 

© 

O 

CM 

CM 

CM 

CM 

© 

© 

tn 

m 

r 

fA 

© 

fA 

*r\ 

m 

•& 

at' 

‘•St 

at 

at 

at 

at 

st 

at 

at 

© 

at  CM  rH  lO  rH 


at  CM  rH  ID  rH 


at  CM  OO  CD  OO 

irl  CM  rH  to 
©  ©  ©  oo 
©  ©  r-  co 

ID  ID  00  OC 

HHKHO 

tO  UD  ID  U3 

in  ©  in  m 
pH  rH  O  O 
©  tn  ©  © 
in  in  m  m 


rt  cm  ©  at  © 
in  in  m  tn  © 


55 


1015  672  872 


TABLE  IV  -  CONTINUED 


[ 


in 

00 

CTl 

10 

Ol 

ID 

UO 

00 

jt 

ID 

3 

in 

CTl 

CO 

00 

ID 

CO. 

m 

in 

00 

CTl 

fs» 

ID 

CTl 

in 

in 

3 

Ki 

ID 

fA 

fH 

cn 

rH 

fA 

K\ 

o 

3 

ID 

-3- 

rH 

O 

Ol 

o 

K> 

CN 

in 

Ol 

00 

O 

ID 

ID 

O' 

CT) 

cn 

00 

ID 

X 

O 

ID 

f*s 

cn 

O 

o 

3 

rH 

3 

CN 

o 

O 

00 

O 

fA 

00 

o 

ID 

fH 

CN 

3 

O 

o 

o 

o 

vD 

CTl 

UD 

o 

ID 

rH 

iH 

fA 

in 

O 

ID 

m 

rH 

in 

rH 

CN 

ID 

00 

in 

fA 

CN 

ID 

fA 

00 

00 

00 

& 

UJ 

• 

« 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

’  • 

• 

_  * 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

* 

« 

O 

O 

rH 

o 

rH 

rH 

rH 

rH 

|H 

rH 

cs 

rH 

rH 

CTl 

o 

rH 

fH 

rH 

in 

o 

CN 

o 

O 

CN 

rH 

tH 

rH 

cs 

cs 

o 

fA 

in 

00 

cn 

ID 

CTl 

ID 

ID 

oo 

3 

o 

ID 

•3- 

in 

CTl 

00 

00 

ID 

00 

in 

in 

00 

CTl 

fs. 

ID 

cn 

00 

in 

in 

3 

fA 

ID 

(A 

rH 

cn 

iH 

(A 

fA 

o 

■3 

o 

ID 

3 

tH 

O 

cn 

o 

fA 

CS 

in 

CTl 

00 

A> 

O 

ID 

ID 

o 

Ol 

Ol 

oo 

ID 

o 

O 

ID 

K 

cn 

n. 

O 

O 

3 

tH 

o 

3 

CN 

o 

o 

OO 

•3- 

O 

fA 

fA 

O 

3 

ID 

rH 

3 

CN 

o 

o 

fA 

O 

oo 

rH 

rH 

o 

o 

o 

rH 

rH 

o 

O 

-a* 

rH 

rH 

o 

CN 

o 

O 

rH 

C, 

00 

O 

fH 

fA 

o 

rH 

o 

o 

o 

o 

CN 

O 

o 

O 

o 

CN 

o 

o 

o 

e> 

o 

o 

o 

O 

o 

rH 

o 

o 

o 

O 

o 

O 

O 

O 

o 

O 

o 

cJ 

o 

o 

O 

O 

CC 

O 

rH 

CM 

rH 

CN 

rH 

fH 

rH 

f- 

fH 

CN 

CN 

CN 

rH 

CN 

fH 

rH 

fH 

m 

O 

O 

CN 

CN 

rH 

tH 

fA 

fA 

fH 

3 

* 

1 

rH 

1 

rH 

i 

a 

rH 

rH 

rH 

fA 

rH 

rH 

rH 

»H 

CN 

rH 

rH 

CN 

cy 

0 

CN 

rH 

rH 

rH 

in 

tA 

rH 

rH 

fH 

fH 

cs 

rH 

rH 

tH 

fH 

fH 

CN 

o 

o 

o 

o 

o 

O 

o 

O 

O 

O 

o 

o 

o 

O 

o 

O 

o 

o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

O 

O 

o 

z> 

O 

c 

O 

o 

o 

c 

O 

O 

o 

o 

o 

o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

O 

O 

o 

o 

O 

o 

O 

O 

o 

h 

rH 

rH 

rH 

rH 

rH 

fH 

rH 

rH 

fH 

rH 

tH 

rH 

rH 

rH 

rH 

fH 

rH 

iH 

rH 

rH 

rH 

tH 

fH 

fH 

rH 

fH 

rH 

rH 

rH 

rH 

o 

i 

CC 


H 

< 

X 

o 

X 

DC 


oooooooooooooooooooooooooooooc 

OOOOOOOOOOOOOOOOOOIAOOOOOOOOOOO 

WN40  JOOSO  3  0  aiOIDints.OOOrHCSlDOOOOIDf''CHCH3f'. 

OHHKlHOOONONrllfl  H-rt  OO^MHOOONOOH  HH  tH 

9—4 


ooocoooooooooooooooooooooooooo 

OOOOOOOOOOOOOOOOOOOOOOOOCDOOOOO 
O  3  00  OMOOJOMOOlOMlfltfOlOHrfNOIOONJf^tfifJ 
HNNHCMfMNHNONHCNOHHNOtNNHNHfMHrtNNHKI 
OOOOOOOOrHOOOrHOOOOOOOOOOOOOOOOO 

oooooooooooooooooooo oo oooooooo 


6  00  fs  (M  K\01lfl0)0U0rf(0NONNlCOK,lN00Nl01ACnC0K1C0N0) 
MNJOHKlKlJ  N'lfl  WOHHO«)K\MOHHjtONONHOHN 

U  . . .  i  . . .  • 

ONlOONiinOOUlOQOOOHN  fAO  0fAf~.00000000 
rH  r-H 


r 


c  CC  HHHlDHHHHNHNNrljfNHHHM  rt'H  HHNHrlrlHNH 


i' 

4 


CO  CD  00  SO  00  ID  00  CO  PH  CO  CO  CO  CO  CO  CD  CO  CO  O  CD 

O  CO  00  O  O  COO  CO  O  CO  CO  CO  CO  CO  CO  Mtfi  HU 

X  JtNJtrf  jf  oj  eiN  ounoicnmoi  cn  ooo 

O  lAsMnWNKUflOsftDOOOOOOJ  J-J-  J  JfO  O  OO  IA  ONW 

UJ  •  •«..•••.. . . . .  •  »  .  • 

HHOHOC'JHOOO  OO  OOOOOGOOOOOOOQOOOCX 


OOCflOOOOCOCDCOOCDKlOOOOOCOCOCDCOCDCOOOOOOOCOOOCO 
OCOOOOOCOOOCOOOOOOOlOCOCOCOCDCOOOOOCOCOOiHCO 
O  ^^JJ.tOJOfflSOOOOOClfflOlOlClBlOOOOd'mOOOO 
00  O  OH  OOOONHOJtJJtfaHHHrtrtHNtf^NHHICNO 

UJ  .......  ........ . . 

OOOOOOOOOOOOO  HO ooooooooooooooo 


OHONHHOOHHHHHOOOOOOHHHHOOi 
t  I  I  I.  I  I  I  I  I  I 


HHHHHHHHHHrtHHHHHHHH 


I  H  rH  rH  rH  ?H  i—C  H  i 


OOOOOOOOOOOOO O  OOOOOOOOOOOOOOO 
r>  O OOOOOOOO  0-0  OC  C.OOOOOOOOOOOOOOOO 
<C  i.  •  ••»•«».»•«■<.••»••  ........... 

(—  HHHHrlHHrlHHH  rH  (H  HHHHHHHHHHHHHHHHH 


OOOOOOOO  OO  OOOOOOOOOOOOOOOOOOOO 
I  O  OOOOOOOOOOO OOOOOOOOOOOOOOOOOOO 

"  ncoconki  P-.r^ 

c  c  . . 

UJ  OOOOOOOO  0  0.0  OOrHOOOOOOOOOOOOOOtHO 

3 
2 

H  H*  OOOOOOOOOOOOOOOOO  OOOOOOOO  o  ooo  o 

Z  <  OOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

O  X  dT'JCM^-cojs-j-CM^eoM  co, oo  cfi-a- 

CJ  O  HHHHOHHOHOOOOOOHHHHHHOOOHH  HH  r-4  rH 

X  OOOOOOOOOOO  0.0  OOOOOOOOOOOOOOOOO 

I  K  . . . . . . 

OOOOOOOOOOO  0.0  OOOOOOOOOOOOOOOOO 

> 


UJ  N  CD  05  N  C  CO  fOID  ID  O  Lfl  (D  H  CChCD  OOMOOCNCIlCMNOlNOlNmOUUAi 

-4  lACnOlOOHHN^fACflOitOlANNCHAOHj-NHtCOOinHOCMJtO 

CO  u  . . .  .  . 

<  OOJ^r-IOOr-tOJtOfO  POPO  NOSOlJtlAIDIDHNtMlOKMOKlOO 

H  r-(  rH  rH  Kl  CM  rH 


CC  HHHHHHHNHHHHHNNHHHHHHNHHHHHHNH 


<  HHHHHHHWHHHHHNNHHHHHHNHHHHHHNH 


NOHNlAt0NC0O|C|JfU5C0NC0NriOlr4NK5jCNO)KlC8d'0001C5 

CTUAtASNNNMJICJlCHIJ101Cl(nOOOHHHHHHNOCONN« 

cococor^H-i^.H>.r>r^i^r^r^r->!^.r^i^.D-.h»Hs^r^fvrv«t^  p-cd  mrtouo 

2  COCOCOOOOOOOOOOOOOOOOOOOOOOOrHOOCOWOP- 

u.  p*.  p^  p*  p*  p»  a,  p^  a.  a*  p*.  r*.  p«.  p»  p- p>.  r*.  p*  a*  r-.  r^.  ps  ps.  t*»  p^  p*.  r*  co  oo  oo 

cncnocncncncncnmocncncncnocncncncricnincntncncn  cocao  u\  o 
HOHHHHHOHHHHrlHHHHHHHHOHHrlOHJfHjf 
PAOfAOOOOPAOPAOOOOPAO  OO  OOOOOOOOOK5PACA 
IflHlOHHHHUlHl/lHHHHlOHHHHrlHHHHHHHininiA 


”0  lO  N  00  Cl  O  H  CVK\  JUMONCOCOrtMIA  Jtin  CC  P->  CO  C  O  H.PI  (O  Jf  IP 
00  00  CO  OO  C  C  C  CC  cnOIOIOlOlOOOOOOOOOOHHnHHH 

HHHHHHHHHHHHHHrIH 


57 


TABLE  IV  -  CONTINUED 


r> 

in 

O)  10  CM 

in  d 

CT1 

10  10  10 

cm  10  oo  in  m 

10  lO  KI  CM  CM 

(■*>. 

10 

ri  N  00 

id  oo 

M 

OIOU0 

00  10  o  o  o 

10 

ID  N  oo  CD 

X 

CM 

o 

O  M  10 

o  ki 

t*- 

oi  oi  oi 

M  O  d  O  M 

o 

OIKIHO 

O  MOMOMMMOMOOOlOOOOOOdOCON'lN'lMlOOddK'llOf" 

111  . . . . . 

iO0inO(Q0OOin0OHOOOOOOONNHIsOOOO4Nlfl 


NoiAOQiDMomii'omoooisioioocvio  co. in  m  o  ioaio  ki  cm  cm 

NOtOOrlNOOOtOCOOHOOOOIOllfiOOQIDObOO  10  UO  Is"  00  Ol 
O  NOOOOHIDOOMOSOOOOICIOIOHO^OSOWOIWHO 
CO  M4040010l40N4  0N4  40  0H4000004rlH|c>00 

UJ  . . . 

MOOOOMNIOOOOOOOO  CM, CM  OOOOOOOOOO  0.0  O 


CC  HHIOHSHHHIDHHNHHHHHOrKON  r-100  OHOOMMN 
II  I  III  I  I  I  1  I  I  I  I 


O'  WHlAHKUCirvHKAHHHHrHHfOHAHHHHHK\riHHHOH4 


O  O  O  O  0-0  OOOOOOOO  0,0  O  O  O'O  OOOOOOOOOO 
X  0000000.0  0000000000000000000000 

<  . . . 

H  H'H  rHHHHHHHHHHHHHHH 


oooooooooooooooooooooooooooooo 

I  o  oooooooooooooooooo  O'O  OOOOOOOOOO 
X  t'*daidoieooodcidddt''ddoot''ddr''t^o>ddi,'«f*»odc0 
cc  . . , . . . 

NONONNICONriOHO  0.0  WK\OOHOONOOOOOlHK\ 


-H  OOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
<  OOOOOOOOOOOOOOOOO  OOOOOOOO o  o  o  o  o 
X  oneocoooooh't^oooooo  oo  co  Njttfmo^Mcoj’d'cocoood'JtoMco 
O  0.0  in  O  in  o  rH  O  in  NO  NOOOHHHONHrtinOOHrUON^ 
X  OOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

CC  . . . . .  • 

O  O  O  OO  00000000000000000000.0  0000 


HHOMDHH40NNONI»COiniOb0444NNin«)01NU)NN 

HCO4'inHinHCO|CiC1HtOf''J'KAClNHSSHinONNtn,0]NNO 

U  . . . . . . .  . 

0N0KA04  cki  or^j-csiooooomooooooencMCMrHoo 

iH  CM  Cl 


GO  U3HHHHOOOOHHHHHNNN  d  ID  rlHHHHHrlrlHHKIHN 


<  ID  H  l — IHHcOOOHMrHHHNNC'JJ'UjHHHHHHHHHHt^HN 


in  Cl  Ol'A  iniDCCOlHIDONHNKMDMniDSOOOIONinNt^eOlDH 
cm  o  mcm  CMnJCMCMfniniDiDr^r^n-oooooooooooooooiCTiCTi.i-J-oijD^' 
KI  O  OO  OOOOOOOOOOOOOHHHHHHHHNNICO14 
NKUflKlKUflWWIAKMfllOKMnm  KAWV1  K\MK\KMOIAK\KMOIfllOd' 

Z  0)00000000000000000000000  OO  O  CM  d  d 
CO  KlddddddddddddddddddddddddddddUl 
ll  oooooocooooooooococooooooooocooooooooooocooooocacoooooooeoco 

in  in  o  in  o  in  tnin  oinintntnininminintnintninooinininininin 
rlrlKArl  KIM  MMKIr*IMMMMMMMMMMMrHKlMMMMMM'r"l 
K)0K)0K)000K1000KAK)K)0000CC)K)K\K)0C0;Mh)O 
if,  M  inHinHHrlinHHHinmtAHHHrtHHlDlMnHHHiniilH 


“3  1DNMOOHNM  4in  ID  f''00  O)  O  M  IN  KV.it  irilONOOOlOHN  Kid  in 
M  M  M  M  CM  CM  CM  CM  CM  CM  CM  CM  CMCMKIKIKIK’linKIK'IKIKINIdddddd 
M  M  M  MM  MMMMMMMMMMM  MrH  MMMMMMMMMMMM 


58 


9 


CM 

o 

to 

O 

CM 

Ol 

to 

to 

to 

o 

Jj- 

CM 

Cl 

00 

Cl 

00 

to 

CM 

rH 

fA 

cn 

on 

Cl 

rH 

00 

C* 

tn 

CO 

X 

O 

o 

o 

a» 

a» 

O 

<T> 

cn 

Cl 

00 

m 

00 

o 

st 

o 

a* 

o 

fA 

CM 

CTi 

to 

o 

o 

o 

00 

o 

to 

a~ 

UJ 

• 

• 

• 

♦ 

• 

• 

• 

• 

* 

* 

• 

-  • 

• 

• 

m 

in 

CM 

tA 

rH 

tA 

O 

o 

o 

o 

o 

o 

fA 

o 

i 


C  OIOONOHDIDIDOJ-NOIOQ 
aiOOlOCMrHfACncnairHeo.3-fAlO 
O  ooocMf^oaiaiCTiootAcoo-a- 
00  O  O  O  rH  A>  OO  OONNNON 

UJ  . . . 

OOOONONNNCONOO 


cc 

7 

6 

2 

5 

6 

H 

1 

I 

I 

I- 

I- 

“1 

4 

1 

a 

4 

3 

1 

4 

5 

rH  tn 

3 

3 

i 

1 

3 

1 

3 

rH 

o  o  o  o  o 

o  o 

o  o  c  o 

o  o  o 

=3 

to  o  o  o  o 

o  o 

O  CD  O  O 

o  o  o 

< 

•  • 

•  •  •  • 

•  •  • 

H 

1 

1 

1 

1 

1 

rH  rH 

I 

t 

T 

T 

I 

I 

I 

o  o  o  o  o 

o  o 

o  o  o  o  o  o  o 

1 

o 

o  o  o  tn  o 

o  o 

o  o  c  o 

o  o  o 

X 

00OMAW 

•=f,  o 

OOU1J 

M  jy  O 

a 

cc 

•  • 

•  •  •  • 

•  •  • 

UJ 

3 

3 
0 

4 
4 

rH  tn 

3 

3 

1 

1 

3 

1 

3 

Z3 

rH 

. 

Z 

— 

h 

H 

O  O  OO  O 

o  o 

O  O  O  CD 

o  o  o 

z 

< 

O  O  O  O  O 

o  o 

O  C  OO 

o  o  o 

o 

X 

ID  O  J'Ol  S 

00  o 

OOW5 

co  oo  o 

a 

o 

A*  ID  rH  CM  fA  CM  rH 

H  H  H  rH 

OWN 

X 

O  O  O  O  O 

o  c 

o  o  o  o 

O  O  O 

• 

a: 

•  « 

*  •  •  • 

•  •  * 

O  O  O  O  O  O  O 

o  o  o  o 

0 

0 

0 

> 

UJ 

lAONKl  J 

-3“  -=r 

•3-  -3-  A»  IA 

moo 

_j 

IA  O  O  O 

o  o 

o  o  o  o 

tn  cm  tn 

ca 

a 

•  s  * 

•  *  4  • 

•  •  • 

< 

o  o  o  o  o 

o  o 

o  o  o  o 

o  o  o 

00 

rH 

rH 

rH  tA 

oo 

rH 

CD 

CD 

CO 

CM 

CM 

00  rH 

tn 

< 

rH 

iH 

rH  fA 

00 

rH 

CO 

CO 

CO 

CM 

CM 

00  rH' 

tn 

tn 

cn 

to  oo 

CM 

St 

o 

rH 

CM 

tn 

00 

tn  in 

.st 

CM 

CO 

cn  K 

to 

to 

r^. 

r- 

r-. 

St  to  A* 

St 

oo 

r> 

a»  tA 

to 

to 

CO 

CO 

co 

CO 

^ t 

tn  r*> 

St 

o 

cn 

cn  cn 

to 

to 

CO 

CO 

CO 

CO 

h-  00 

to 

z 

CO 

in 

oo  oo 

fA 

fA 

tn 

tn 

tn 

m 

tn  oo 

O 

00 

cn 

^ t 

•a- 

00 

00 

00 

oo 

00 

00 

00 

00  oo 

CT1 

u. 

CO 

cn 

cn  cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn  cn 

Cl 

in 

in 

o  m 

tA 

IA 

in 

tn 

tn 

in 

in 

in  in 

m 

rH 

rH 

rH  O 

o 

o 

o 

o 

o 

o 

o 

OH 

o 

tn 

o 

iH  tA 

fA 

fA 

tn 

tn 

tn 

m 

tn 

tn  tn 

fA 

in 

rH 

fA  m 

IA 

in 

in 

tn 

in 

tn 

in 

in  in 

IA 

”3 

to 

A. 

oo  a> 

O 

rH 

CM 

fA 

St 

in 

ID  A-  00 

Cl 

.3- 

St  St 

in 

tn 

in 

in 

in 

tn 

in 

tn  in 

in 

rH 

rH 

rH  rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH  rH 

rH 

59 


This  (Q,r)  implies  that  the  part  is  not  stocked  by  the 
maintenance  unit  but  is  ordered  from  the  source  of  supply 
as  demands  occur. 

The  data  show  that  expected  costs  of  on  hand  inventory 

rise  very  steeply  as  the  required  level  of  A  approaches 

s 

unity. 


TABLE  V 


EXPECTED 

COST  OF  ON  HAND 

INVENTORY  VERSUS 

REQUIRED  A 

s 

A 

s 

E(C) 

No.  of 
Iterations 

No.  of 
Negative 
Safety 
Levels 

0.90 

$16.80 

31 

156 

0.91 

$22.02 

64 

151 

0.92 

$38.14 

67 

143 

0.93 

$79.57 

140 

123 

0.94 

$168.64 

141 

104 

0.95 

$418.04 

295 

73 

0.96 

$645.76 

238 

41 

0.97 

$1225.36 

359 

19 

0.98 

$1619.00 

322 

7  ' 

0.99 

$2406.83 

386 

2 

Figure  1  plots  the  expected  cost  of  on  hand  inventory 

E  (C)  as  a  function  of  A  . 

.  s 
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Figure  1.  Expected  Cost  of  on  Hand  Inventory  E(C) 
as  a  Function  of  Required  Expected  Fraction  of 
Unit  Systems  Up 


2.  Expected  Cost  of  On  Hand  Inventory  Versus  Required 
Number  of  Unit  Systems  Up  with  an  Assurance  of  90 
Per  Cent^ 

Program  PKSMIN  was  run  for  six  values  of  k  =  45, 
...,50  at  a  required  assurance  of  90  per  cent  in  order  to 
develop  a  functional  relationship  between  the  required  value 
of  k  and  the  expected  cost  of  on  hand  inventory  at  this 
assurance  level.  Table  VI  summaries  the  results. 

TABLE  VI 


EXPECTED  COST  OF  ON  HAND  INVENTORY  VERSUS  REQUIRED 
NUMBER  UNIT  SYSTEMS  UP  WITH  90  PER  CENT  ASSURANCE 


k 

E(C) 

No.  of 
Iterations 

No.  of 

Negative 

Safety 

Levels 

45 

$29.07 

103 

150 

46 

$84.97 

102 

122 

47 

$387.88 

177 

85 

48 

$1117.47 

218 

49 

49 

$2235.32 

278 

11 

50 

$3519.92 

303 

0 

Figure 

2  plots  the  expected  cost  of  on 

hand  inventory 

E  (C)  as  a 

function  of  the 

required  number 

k  of  unit 

systems  up 

with  90  per  cent 

.  assurance. 
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Figure  2.  Expected  Cost  of  on  Hand  Inventory  E(C) 
as  a  Function  of  the  Required  Number  k  of 
Unit  Systems  Up  with  90  Per  Cent  Assurance 


V.  OBSERVATIONS  AND  CONCLUSIONS 


-  A.  OBSERVATIONS 

1 .  Marginal  Analysis  may  Produce  only  a  Near-Optimal 

Solution 

The  marginal  analysis  procedure  presented  in 
Section  III  sometimes  overshoots  the  optimum  solution. 

That  is,  the  r  vector  found  using  the  procedure  sometimes 
results  in  a  higher  than  required  availability  and  a  higher 
expected  cost  of  on  hand  inventory  than  needed  to  satisfy 
the  constraint.  This  behavior  of  marginal  analysis  was 
observed  during  the  hand  calculation  of  the  solution  to 
a  small  problem  with  four  unit  systems  and  four  part  types 
where  the  constraint  was  E(Z)  £.3.6.  Marginal  analysis 
yielded  an  E(Z)  of  3.709  and  an  E (C)  of  $53.87.  It 
was  possible,  however,  on  the  last  iteration  to  achieve 
an  E (Z )  of  3/601  at  a  cost  of  $48.63  by  increasing  a 
different  r^  than  the  one  associated  with  the  maximum 
Aj.  A  dynamic  programming  procedure  would  overcome  this 
shortcoming,  but  dynamic  programming  rapidly  becomes  compu¬ 
tationally  cumbersome  as  the  number  of  decision  variables 
increases.  In  the  sample  problems  discussed  in  Section  IV 
there  were  159  decision  variables,  and' an  many  possible 
applications  of  tKe  model  the  number  would  be  in  the 
thousands.  Thus  the  simple  marginal  analysis  procedure  may 
be  more  computationally  feasible  than  dynamic  programming 
for  many  applications  of  interest.  Further,  the  relative 
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differences  between  the  optimal  solution  and  the  marginal 
analysis  solution  tends  to  decrease  as  the  size  of  the 
problem  increases. 

2.  Program  PKSMIN  is  Much  Faster  than  Program  EZMIN 
Sample  problems  1  and  2  in  Section  IV  resulted  in 

about  the  same  values  of  E(Z)  and  of  E(C).  For  problem 

1  E (Z)  was  47.58  and  E(C)  wa's  $418.04.  For  problem  2 
E  (Z)  was  47.46  and  E(C)  was  $387.88.  On  the  IBM  360 
program  EZMIN  required  over  four  minutes  to  compute  the 
results  for  problem  1,  while  program  PKSMIN  solved  problem 

2  in  nineteen  seconds.  The  reason  for  this  difference  is 
that  E (Z)  a  more  complex  function  than  Pr(Z  >  k) .  For 
a  given  r  vector  E(Z)  is  the  sum  over  k  from  one  to 
S  of  Pr(Z  >  k) .  A  comparison  of  the  JMAX  segments  of 
subroutines  OPTIMR  and  OPTMZ  shows  that  OPTIMR  requires 
more  computations  at  each  iteration  (See  Flow  charts  B4a 
and  C3a. ) 

3.  Functional  Relationships  of  Expected  Cost  Versus 
Rgyired  Availability  are  Readily  Obtained 

The  functional  relationships  of  expected  cost  of 

f 

on  hand  inventory  versus  availability  discussed  in  Section 
IV  were  obtained  by  simply  adding  a  DO  loop  to  the  two 
main  programs.  Thus  for  E(C)  versus  A  the  value  of 
A  was  initially  set  to  0.89,  and  at  each  iteration  of  the 
DO  loop  Ag  was  increased  by  0.01  for  ten  iterations. 
Similarly  for  E (C)  versus  k  at  the  90  per  cent  assurance 
level  k  was  initially  set  at  44  and  incremented  by  one 
at  each  iterations  for  six  iterations. 
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Using  the  functional  relationships  of  E(C)  versus. 

availability  we  can  estimate  the  maximum  availability  which 

can  be  achieved  if  there  is  a  constraint  on  E(C).  Suppose, 

for  example,  we^  desired  to  maximize  A  subject  to  E (C) 

<  $1000.00.  Prom  Table  V  in  Section  IV  we  see  that  the 

maximum  A  would  be  between  .96  and  .97.  To  determine 
s 

maximum  Ag  more  precisely  and  to  determine  the  associated 
r  vector  we  could  run  program  EZMIN  for  a  series  of  values 
of  A  ,  'say  .961,  .962,  ...,  .969.  The  resulting  E(C) 

5 

for  one  of  these  A  values  should  be  very  close  to 

b 

$1000.00. 

Computing  time  in  developing  the  E (C)  versus  avail¬ 
ability  functions,  can  be  conserved  if  at  each  iteration 
after  the  first  the  starting  values  of  the  r^  are  'the 
final  values  computed  on  the  previous  iteration.  This  can 
be:  accomplished  if  the  initializing  subroutines  INITAL  and 
INTLK  are  skipped  on  the  second  and  subsequent  iterations. 
This  technique  reduces  the  computing  time  to  about  half 
that  needed  if  the  initializing  subroutines  are  called  at 
each  iteration.  ;  . 

»  x 

B.  CONCLUSIONS 

i .  Advantages  of  the  Model 

The  model  provides  the  first  technique  of  which 
the  author  is  aware, 'for  explicitly  taking  into  account 
system  availability  in  the  determination  of  continuous 
review  inventory  policies  for  repair  parts.  The  model 
allows  for  variable  safety  levels  and  yields  a  set  policies 
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where  the  ratio  of  marginal  availability  to  marginal  cost 
of  inventory  investment  is  approximately  the  same  for  all 
part  types.  The  result  is  that  the  model  tends  to  favor 
low-cost,  high-demand  items  with  relatively  high  safety 
levels.  High-cost,  low-demand  items  are  less  well  protected. 
All  items  are,  however,  well  enough  protected  to-  insure, 
under  the  assumptions  of  the  model,  that  the  availability 
requirement  will  be  met. 

Cpmpared  to  models  which  might  .more  realistically 
represent  the  "real  world"  the  one  presented  here  possesses 
the  advantage  of  computational  feasibility.  Programs 
EXMIN  and  PKSMIN  produce  solutions  in  a  reasonable  amount 
of  time.  Refinements  in  the  programming  could  probably 
improve  run  times. 

2 .  Limitations  of  the  Model 

The  limitations  stem  primarily  from  the  assumptions 
which  were  made  for  the  sake  of  mathematical  simplicity. 

For  example,  the  constant  lead  time  assumption  is  not 
usually  true.  Lead  times  can  have  considerable  variability. 
Le.  t  time  variability  will  tend  to  decrease  the  actual 
availabilities  achieved  compared  to  those  predicted  by  the 
model.  On  the  other  hand  the  model  does  not  consider 
the  fact  that  requisitions  for  parts  which  are  causing  a 
unit  system  to  be  down  have  a  higher  -priority  and  thus  a 
shorter  lead  time  than  normal  replenishment  requisitions'; 
This  fact  will  tend  to  make  achieved  availabilities  higher 
than  those  predicted. 
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The  optimal  cannibalization  assumption  is  seldom 
completely  true.  The  significance  of  this  assumption 
decreases,  however,  as  the  required  availability  increases, 
since  at  higher  availability  levels  we  expect  fewer  back- 
orders  and  therefore  fewer  occasions  for  cannibalization. 

The  significance  of  the  assumption  of  zero  or  one  states 
for  unit  systems  will  depend  upon  the  stringency  of  the 
serviceability  criteria  used  to  determine  whether  a  unit 
system  is  up  or  down.  The  more  stringent  the  criteria, 
the  more  parts  will  have  b^  equal  to  a^.  Parts  with 
bj  equal  to  a^  will  generally  be  better  protected  than 
parts  with  b^  less  than  a \ .  The  most  stringent  criteria 

possible  would  be  that  every  part  on  the  unit  system  must 
be  working  if  it  is  to  be  counted  as  being  up.  Thus  we. 
can  see  that  the  inventory  policies  computed  by  this 
model  and  the  resulting  expected  costs  of  on  hand  inventory 
are  highly  dependent  upon  the  serviceability  criteria  tor 
the  unit  system. 

•  The  assumption  that  parts  are  replaced  immediately  with 
zero  lead  timemakes  it  necessary  to  correct  the  avail¬ 
abilities  computed  by  the  model  for  down  time  due  to  finite 
time  to  repair. 

3.  Uses  of. the  Model 

The  model  was  developed  for  use  in  determining 
direct  support  inventory  policies  for  repaire  parts  for  a 
given  population  of  a  particular  major  item. 
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Another  use  which  suggests  itself  is  determining 
optimal  maintenance  floats.  Suppose  it  were  required  that 
at  least  k  unit  systems  be  up  with  probability  .90. 

Assume  that  the  capability  of  the  maintenance  unit  to  make 
repairs,  given  that  needed  parts  are  available* is  fixed. 

Then  the  probability  that  k  unit  systems  are  up  is 
determined  by  the  total  number  S  of  unit  systems  and 
the  inventory  policies  for  repair  parts.  The  quantity 
S  -  k  is  called  the  maintenance  float.  An  interesting 
problem  is  determining  the  optimal  level  of  the  maintenance 
float.  The  model  presented  in  this  thesis  could  be  used  to 
estimate  the  optimal  float  level  as  follows.  Compute  the 
expected  cost  of  on  hand  inventory  needed  to  assure  with 
probability  .90  that  k  unit  systems  are  up  with  zero;  one* 
two,  ...  unit  systems  in  the  maintenance  float.  Then  choose 
that  level  of  maintenance  float  which  minimizes  the  sum  of 
the  cost  of  the  float  plus  the  on  hand  inventory. 

4 .  Extensions 

A  number  of  possible  extensions  of  the  model 
appear  to  be  worthy  of  investigation.  For  example,  relax¬ 
ation  of  the  assumption  that  all  unit  systems  are  identical 
would  be  useful.  A  next  higher  level  of  complexity  of  the 
system  would  be  one  in  which  there  are  two  types  of  unit 
systems.  For  example,  one  type  of  unit  system  might  be  a 
particular  type  of  weapon.  The  other  unit  system  might  be 
a  piece  of  fire  control  equipment  which  controls  the  fires 
of  several  weapons.  If  a  fire  control  unit  is  down  all  the 
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weapons  it  controls  are  also  down.  Application  of  the 
approach  of  this  model  to  optimal  inventory  policies  for 
such  a  system  seems  to  offer  interesting  possibilities. 

In  any  case  it  is  hoped  that  the  model  presented  here 
will  prove  useful  in  the  development  of  system  availability 
oriented  operating  policies  for  military  repair  parts 
inventories . 


70 


APPENDIX  A 


CONVENXITY  CONSIDERATIONS 

In  order  to  insure  the  optimality  of  the  r  vector 
found  using  the  procedure  outlined  in  Section  III  it  is 
necessary  that  the  cost  function  be  convex  and  that  the 
region  in  r-space  defined  by  the  constraint  be  convex. 

1.  The  Cost  Function 

The  objective  function  in  the  model  presented  in 
this  paper  is  the  expected  cost  of  on  hand  inventory, 
which  is  given  by 

N 

D (C)  =  l  C.D  (Q  •  ,r  . )  . 
j=l  J  J  J 

Since  the  Qr.,  are  fixed  by  the  operating  level  of  supply 
policy  as  described  in  Section  III,  we  are  concerned 
only  with  convexity  with  respect  to  the  r j ,  j  =  1,...,N. 
Further,  since  E(C)  is  a  sum  of  terms  each  of  which  is 
a  function  only  of  r  j ,  E(C)  is  convex  if  each 
is  convex.  Since  Qj  is  fixed  let  us  drop  it  as  an  argu¬ 
ment;  also  let  us  fix  j  and  drop  the  subscript.  Hadley 
and  Whitin  [Ref.  8,  p.  184]  show  that 

D(r)  =  (Q  +  l)/2  +  r  +.-  B (Q,r)  , 

where  B (Q,r)  is  the  expected  backorders,  given  a  (Q,r) 
policy.  D(r)  is  linear  in  r  except  for  the  B (Q,r)  term. 
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Thus  D(r)  is  convex  if  B(Q,r)  is  convex  in  r.  Again, 
since  Q  is  fixed,  let  us  drop  it  as  an  argument.  Hadley 
and  Whitin  [Ref.  8,  p.  184]  show  that 


B  (r)  =  (1/Q)  l  y  [P  (y+r+1; y)  -  P  (y+r+Q+1; y)  ],. 
y=0 


Let 


ArB(r)  =  B  (r+1)  -  B{r)  „ 

and 

A^B(r)  =  ArB  (r+l)  -  ArB  (r)  . 

2 

D (r)  is  convex  if  ArB(r)  1  0  for  all  r.  Since  the  par¬ 
ameter  y  is  fixed  for  fixed  j ,  let  us  further  simplify 
the  notation  by  letting 


P(x)  =  P (x; y)  , 


p (x)  =  p (x;y) . 

Further,  let 

B*(r)  =  QB(r)  . 

2 

Then  D(r)  is  convex  for  r  >.  r*  if  ArB*  (r-)  k.  0  for  all 
r  >.  r*.  Now 


ArB*<r) 


r  v  P  (y+r+2)  -  P  (y+r+Q+2) 
y=0  y-P(y+r+l)  +  P  (y+r+Q+1) 


=  2  y[p (y+r+Q+1)  -  p (y+r+1)]. 

y=l 
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Note  that  we  nave  dropped  the  term  for  y  =  0  since  it 
contributes  nothing  to  the  sum. 

Now 


A 


’2 

r 


B*  (r) 


00 


E 

y=l 


p (y+r+Q+2) 
-p(y+r+Q+l) 


-  p (y+r+2) 
+  P  (y+r+1) 


=  £  y [Ap (y+r+Q+1)  -  Ap (y+r+1)] 

y=i 


where 


Ap  (x)  =  p  ( x+ 1 )  -  p  (x) 


u(ittl)  e'11 
(X+1)  ! 


fu  -  x  -  1 
1  x+1 


]  p  (x)  . 


Thus 


A^B* (r) 


Z 

y=i 


y[ 


y-y-r-Q-2 

y+r+Q+2 


P (y+r+Q+1)- 


y-y-r-2 

y+r+2 


p  (y+r+1)  ] 


=  Z 

y=i 


y[ 


y+r+2-y 

y+r+2 


p (y+r+1) 


y+r+Q+2-y 

y+r+Q+2 


p  (y+r+Q+2) 


(r+3-y ) ._ 
r+3-  P(r+2) 


(r+Q+3-y) 

r+Q+3 


p (r+Q+2) 


+  il|±^uip(r„,  -  2M2|iaiP(r+Q+3) 


+  (Q+l)-^||i|ip(r+Qi-2)  -  ... 
+  ... 
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=  V  y(£^|^L)p(y+r+l)  +  Q \  (£^piJi)p(r+y+l). 

Nov;  every  term  in  this  series  is  non-negative  if  r  >  {i— 3 . 
Thus  D(r)  is  convex  in  r  for  r  =  y-3.  Furthermore, 
numerical  calculations  show  that  for  values  of  y  up  to 
40  that  if  Q  >  y*5  then  D(r)  is  convex  for  r>.  -1. 

This  numerical  result  is  probably  due  to  the  fact  that  the 
mode  of  the  Poisson  probability  mass  function  is  approxi¬ 
mately  y  -1,  which  means  that  the  most  negative  values 
of  (r+y+2-u„/ (r+y+2)  in  the  series  are  multiplied  by 
small  values  of  p(x;y).  When  r  -1,  all  terms  in  the 
series  for  y  >.  y  -  1  will  be  positive.  Thus  it  is 
likely  that  the  sum  of  the  positive  terms  in  the  series 
is  greater  than  the  absolute  value  of  the  sum  of  the  terms 
for  which  y  4  1  -  y  is  negative. 

Thus,  as  far  as  the  cost  function  is  concerned,  we  can 
be  sure  that  the  function  is  convex  in  r  for  all 
r  >.  y  -  3,  'where  y  =  y^,..„  ,yN,  and  we  have  reason  to 
believe  the  function  is  convex  for  r  >  -  1  when  Qj^yj-,5, 
j  =  1, . . . ,N. 

2 .  The  Constraint  Function 

The  two  alternate  forms  of  the  constraint  functions 
are,  first,  a  required  expected  number  of  operational  unit 
systems,  expressed  by 

S  N  (Sa-j-kb-i) 

2  n  .2  J  V .  (y.;r)  >  AS  , 
k=l  j=l  y.=0  J  J  s 
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and  second  a  required  probability  bhat  at  least  k  unit 
systems  are  operational,  expressed  by 


N  Sa-s-kb.: 
n  EJ  JY.  (y;r.) 
j=l  y.=0  3  3 


=  P. 


mm 


(k)  . 


If  the  constraint  region  is  to  be  convex,  the  above  func¬ 
tions  must  be  concave  functions  of  r  =  (r.,...frN).  Note 
that  the  left  hand  side  of  the  first  constraint  is'  the  sum 
over  k  of  terms,  of  the  same  form  as  the  left  hand  side 
of  the  second  constraint.  Thus,  if  we  can  show  the  con¬ 
ditions  for  concavity  of  the  left  hand  side,  of  the  second 
constraint  for  a  general  k,  we  can  easily  show  the 
conditions  for  concavity  of  the  first  constraint. 

The  difficulty  in  showing  the  concavity  of  the  second 
form  of  constraint  is  that  it  is  not  a  separable  function 
of  the  r j .  It  can  be  transrormed  into  a  separable  function 
by  taking  the  logarithm,  in  which  case  it  becomes  a  sum  of 
terms,  each  of  which  is  a  function  of  only  one  r ^ ,  i.e., 

N 

In [Pr {Z  >  k}]  =  E  ln[Pr{Z.  £  k}], 

j  =  l  3 


where 


Let 


N 

f  (r)  =  Pr{  Z  >  k? r} ,  =  II  f .  (r. )  , 

j=i 3  1 

f j (rj)  =  Pr{Z j  >  k;r j } , 

4 

I 

<p  (r)  =  In  f  (r)  , 

W  =  ln  fj(rj>- 

Now  it  can  easily  be  shown  that  if  fj(r^)  is  concave, 
<j>.(r.)  is  also  concave,  and  thus  <j>  (r)  is  concave.  Un- 

•J  J 

fortunately,  the  concavity  of  <j>  (r)  does  not  imply  the 
concavity  of  f (r) .  Now  consider 

Sa-i-kb^ 

f,(r.)  =  J4' .  (y  •  ;r . ) . 

J  J  y^=0  J  J  J 

Let 

M  =  Sa^  -  kb j , 

then 

f.  (r.)  =  ¥. (0;r.)  +  ...  +  Y.(M;r.). 

3  J  3  3  3  3 

In  section  III  we  showed  that 

rj+l 

V0;r.)  =  l 

J  J  y.=0  J  J 

*3 

^ j  Cy j ;  )  =  Vj(yj  +  +  1;  -  1)  ,  y^  >  0. 
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Thus  we  have 


f 


j 


rj+1 

Z  ¥ .  (u;-l) 
u=0  J 


M 

+  E  ¥.(y.+r.+l;  - 1) 
y=l  j  3  3 


M+r^+1 

=  £  Y^(u;-1). 

u=0  3 


Now  let 


Vs(V =  W11  -  W' 


then 


M+r-+2  M+r.+l 

Ar  f(r.)  =  S  J  Y.(u;-1)  -  S  3  4'. (u;-l) 


u=0 


u=0 


Let 


=  4,j (M+r..+2;-l) , 


Ar.fi  =  Ar  f (r-j+1>  -  A  f (r.) 
3  j  j  J 


then 


Ar  fj  =  4* .  (M+r+3;-l)  -  y  .  (M+r+2;-l) 
i  J  J  J 


=  (1/Q), 


P  (M+r  j+3) 
-P (M+r+2) 


-  P (M+Qj+r j+3) 

-  P(M+Q+r+2) 


=  (1/Q) [p (M+Q+r+2)  -p (M+r+2)] 


Now  f j  is  concave  for  all  rj  >.  r*  if  the  expression 

in  brackets  above  is  <  0  for  all  r..  >  rj.  Thus  we  seek 

rj  such  that  for  all  r .  >  r* 

3  3“D 
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p (M+r ^+2)  >  p(M+Qj+rj+2) 


or 


y 


1 


(M+r .+2) 

e 

(M+r ^+2)  ! 


(M+Q.+r.+2) 

-L _ _ _ 

(M+Q  .+r  .+2) 
D  D 


or 


Qn  (M+r .+Q .+2) ! 
u  J  <  J  J 

(M+r j+2 )  !  * 

The  right  hand  side  of  -the  above  inequality  is  an  increasing 
function  of  r j ,  Q.. ,  and  M.  It  is  usually  true  that  the 

operating  level  _  olicy  used  within  the  Army  is  at  least  as 
many  days  of  supply  as  the  procurement  lead  time.  This 
means  that  >.  u...  Assume  this  is  true.  Now  if  k  =  S 

and  a.  =  b . ,  then  M  =  0.  If  f.  (r. ;k=S)  is  concave, 

•J  J  J  J 

fj(*j;k<S)  will  also  be  concave,  since  if  the  inequality 
holds  for  M  =  0,  it  will  hold  for  M>0.  Let  M  =  0.  Then 
the  inequality  becomes 


(r  .+Q.+2)  ! 

<  - J - J _ 

(r.+2)! 


For  example,  if  r.  =  -  1  and  Q.  «  1,  then  \i .  must  be 

j  j  j 

equal  to  or  iess  than  2.  In  other  words,  for  u^  <  2,  f  j  (rj) 
is  concave  in  r^  for  all  r^  >  -  1,  Qj  >_  1,  arid  k  <L  S. 

To  insure  the  concavity  of  fj(r^)  for  y  >2,  in  the 
FORTRAN  program  for  the  E(Z)  =  AS  form  of  the  constraint, 

b 
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the  subroutine  INITAL  finds  an  initial  set  of  such  that 

the  concavity  conditions  for  k  =  S  are  satisfied  for  r ^ 

equal  ito  or  greater  than  the  initial  r^.  Similarly,  in 

the  FORTRAN  program  for  the  Pr(Z=k)  =  form  of  the 

constraint,  subroutine  INTLZ  finds  an  initial  set  of  r^ 

such  that  each  Pr(Z_.  =  k;  r^)  is  a  concave  function  of 

r j  for  r^  equal  to  or  greater  than  the  initial  value. 

The  solution  procedure  employed  in  the  computer  programs 

thus  insures  that  each  f^(rj)  is  concave  in  r j  in  the 

region  that  the  marginal  analysis  procedure  searches  for 

N 

an  optimum.  The  question  of  the  concavity  of  f(r)=  II  f.(r.) 

“  j=l  J  J 

remains  open.  To  date  the  numerical  results  have  not  inid- 
cated  that  there  is  a  problem,  i.e.,  the  numerical  results, 
are  quite  reasonable.  No  proof,  however,  of  the  concavity 
of  f  (r)  has  been  found. 


APPENDIX  B 


PLOW  CHARTS  FOR  FORTRAN  PROGRAM  EZMIN 

FORTRAN  program  EZMIN  was  developed  to  solve  the  problem 
of  finding  the  vector  r  which  minimizes  the  expected  cost 
of  on  hand  inventory  subject  to  a  required .expected  number 
of  unit  systems  up.  This  appendix  contains  flow  charts  for 
the  main  program  and  the  following  subprograms: 

Subroutine  PSITAB,  which  computes  a  table  of  values  of 
PI ( Jf  Y)  for  j  =  1 , . . . ,N ,  and  Y  =  1,...,100,  where 

PI  (J ,  Y)  =  ^(Y-ljr.j  =  -  1)  . 

Subroutine  IHITAL,  which  computes  a  set  of  initial 
values  of  r  ^ ,  such  that  E(Z^)  >.  j  =  1,..,,N, 

Subroutine  OPTIMR,  which  uses  a  marginal  analysis  tech¬ 
nique  to  find  an  optimal  r  vector. 

Function  PRZJK,  which  computes  Pr{Zj  =  k;  r^}. 

Function  EBO,  which  computes  ECY^;  rj)* 

Function  PPT,  which  computes  individual  terms  of  the 
Poisson  probability  distribution. 

Function  PCUMT,  which  computes  complementary  cumulative 
terms  of  the  Poisson  probability  distribution. 


Figure  Bl.  Flow  Chart  for  EZMIN  Main  Program 
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COMPUTES  INITIAL  r  .  ,  j  =  SUCH  THAT 

E,(Z.)  =  AZ*S  AND  Pr(Z^3=  S;r.:)  IS  CONCAVE  IN  r. 
D  J  J  j 


COMPUTES  r  VECTOR  SUCH  THAT 
IS  MINIMIZED,  SUBJECT  TO  E(Z)  =  AZ*S 


Figure  B2.  Flow  Chart  for  Subroutine  PSITAB 
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Figure  B2.  Continued 
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Figure  B3.  Flow  Chart  for  Subroutine  INITAL 
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Figure  B4.  Flow  Charb  for  Subroutine  OPTIMR 


Figure  B4 .  Continued 
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I 


Figure  B4a.  Flow  Chart  for  JMAX  Segment  of  OPTIMR 
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MU  =  RHO*T 


Figure  B7.  Flow  Chart  for  Function  PPT 
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Figure  B8.  Flow  Chart  for  Function  PCUMT 
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APPENDIX  C. 


PLOW  CHARTS  FOR  FORTRAN  PROGRAM  PKSMIN 

FORTRAN  program  PKSMIN  was  developed  to  solve  the  prob¬ 
lem  of  finding  the  r  vector  which  minimizes  the  expected 
cost  of  on  hand  inventory  subject  to  a  required  probability 
that  at  least  k  unit  systems  are  up.  This  appendix  con¬ 
tains  flow  charts  for  the  main  program  and  the  following 
subroutines : 

Subroutine  INTLZ ,  which  compute  a  set  of  -initial  values 
of  r^  such  that  Pr{Zj  >.  k)  >.  P(Z,  >.  k)m^n,  .  ,N. 

Subroutine  OPTMZ,  which  uses  a  marginal  analysis' tech¬ 
nique  to  find  an  optimal  r  vector. 

S 

Program  PKSMIN  also  uses  subroutine  PSITAB  and  function 
subprograms  PRZJK,  EBO,  PPT,  and  PCUMT,  flow  charts  for 
which  are  in  Appendix  B. 
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READ:  S  =  NO.  OF  UNIT  SYSTEMS 

KS  =  REQUIRED  NO.  OPERATIONAL 
UNIT  SYSTEMS  > 

N  =  NO.  OF  PART  TYPES  / 

\  PKSMIN  =  REQUIRED  ASSURANCE/ 

\  OL  =  OPERATING  LEVEL  / 

\  ALEVEL  =  ACTIVITY  LEVEL  / 


Figure  Cl.  Flow  Chart  for  PKSMIN  Main  Program 
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INTLZ 


COMPUTES  INITIAL  r^,  j  = 

SUCH  THAT  Pr  (Z  j^_k:  r  j )  >  PKSMIN 
AND  Pr  (Z;  >_KS;rj)  IS  CONCAVE  IN 


OPTMZ 


4 

V 


COMPUTES  r  VECTOR  SUCH  THAT  E (C)  IS 
MINIMIZED  SUBJECT  TO  Pr(Z>K)  >  PKSMIN 


STOP 


) 


A 


A 


Figure  Cl.  Continued 
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N 

COMPUTE  P  (Z>KS)  =  II  Pr(Z.>KS) 

j=i  3 


i 


COMPUTE  EB.  =  EBO  (r . )  =  E  yH*.  (y;r.)  j  = 

J  J  y=l  J  J 


t 


COMPUTE  DTCOST ;  =  C . (1  +  EBO(r.+l)  -  EB.),j  = 

J  J  J 


* 
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Figure  C3.  Continued 
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yes 


* 


PKSTRY  =  Pr (Z>KS) 

PKSTRY  =  Pr  (Z>KS)  [Pr  (Zj>KS)+PI  (J,M)  ]/Pir  (Zj>KS) 


Figure  C3a.  Flow  Chart  for  JAMX  Segment  of  OPTMZ 
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